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ABSTRACT 

We present observations with VLT and HST of the broad emission hnes from the inner 
ejecta and reverse shock of SN 1987A from 1999 Feb. until 2012 Jan. (days 4381 - 9100 after 
explosion). We detect broad lines from Ha, H/3, Mg I], Na I, [O I], [Ca II] and a feature at ~ 
9220 A. We identify the latter Une with Mg II AA 9218, 9244, which is most likely pumped by 
Lya fluorescence. Ha, and H/3 both have a centrally peaked component, extending to ^ 4500 
kms""'^ and a very broad component extending to > 11,000 kms~^, while the other lines have 
only the central component. The low velocity component comes from unshocked ejecta, heated 
mainly by X-rays from the circumstellar environment, whereas the broad component comes from 
faster ejecta passing through the reverse shock, created by the colhsion with the circumstellar 
ring. The flux in Ha from the reverse shock has increased by a factor of 4 — 6 from 2000 to 
2007. After that there is a tendency of flattening of the light curve, similar to what may be 
seen in soft X-rays and in the optical lines from the shocked ring. The core component seen in 
Ha, [Ca II] and Mg II has experienced a similar increase, which is consistent with that found 
from HST photometry. The ring-like morphology of the ejecta is explained as a result of the 
X-ray illumination, depositing energy outside of the core of the ejecta. The energy deposition of 
the external X-rays is calculated using explosion models for SN 1987A and we predict that the 
outer parts of the unshocked ejecta will continue to brighten because of this. We finally discuss 
evidence for dust in the ejecta from line asymmetries. 
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1. Introduction 

During the first decade after the explosion the 
spectrum of SN 1987A was dominated by lines 
from newly synthesized metals in the inner ejecta, 
powered by the decay of ^^Co, ^^Co and "''^Ti. 
Now, more than twenty years after explosion the 
ejecta are involved in an increasingly intense col- 
lision with the circumstellar ring. This is seen 
in both optic al/UV, radio and X-rays (see e.g., 
McCrav l2007l for a review). In the optical the 
collision manifests itself most clearly as a forest of 
increasingly br ight emission lines with velocities of 
~ 300kms~^ ( Pun et al.ll2002 : Groningsson et al 



2008bl ). However, a number of very broad lines, in 
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particular Lya and Ha with velocities of > 10^ 
kms~^, are also evident in the spectrum. Broad 
Ha and Ly a lines from the r everse shock were 
first seen by ISonneborn et al. (Il998l) . It is likely 
that most of that emission is coming from col- 
lisional excitation of the neutral expanding H I 
by shocked electrons, giving an e x tremely broad 
component (' Michael et"aDll998a[ ). iMichael et"ar 



([.1998b, 2003) found from modeling of narrow slit 
observations with HST that this emission was con- 
centrated to a region of ~ ±30° from the equa- 
torial ring, with much less emission from higher 
latitudes. 

The evolu tion of the broad Ha was discussed by 
Smith et al. I ([2005) from a combination of observa- 
tions with STIS and the Magellan 6.5 m telescope. 
An interesting prediction was that the ionizing 
photons from the reverse shock would pre-ionize 
the high velocity ejecta gas and thereby quench 
the broad Lya and Ha emission from the reverse 
shock. Because of the limited signal to noise (S/N) 
and spectral resolution in the STIS spectra and 
differences in the slit width, they could, however, 
not make any firm statement on the evolution of 
the Ha flux. 

Both the Lya and H a evolution f rom 1 999 to 
2004 were discussed by iHeng et al. (l2006l) from 
STIS observations, who found an increase in both 
the Lya and Ha fluxes by factors 5.7 to 9.4 and 2 
to 3, respectively. The larger increase of Lya was 
explained as a result of resonance back-scattering 
of the Lya photons. They also discussed what 
they called the 'interior emission', i.e., line emis- 
sion. This could either originate in low velocity 
ejecta or from charge transfer in the post-shock gas 
of slow H II with high vel ocity H I. A carefu l analy- 
sis of this mechanism by iHeng fc McCrav ( 2007 ) , 
however, found that this 'interior emission', cor- 
responding to the 'broad component' in Balmer 
dominated shocks, was too strong to be explained 
by the charge transfer model. 

[France et al. I (l2010h discussed the Lya and Ha 
emission based on STIS observations from 2010. 
They find a further increase in the Lya and Ha 
fluxes, for the latter by a factor of ~ 1.7 from 2004 
to 2010. Most interesting, from the large extent 
of the blue wing of Lya they propose that Lya 
photons from the hot-spots are boosted in energy 
by scattering against the expanding ejecta. This 
works for Lya since it is a resonance line but not 



for Ha. 

Recently, iLarsson et al. ( 2011 , in the following 
Lll) found from photometry of the ejecta images 
from HST that while the flux from the ejecta de- 
cayed as expected from the ''^Ti decay until 2001 
(day ^ 5000), it has increased since then. By 2010 
the flux was ^ 3 times higher in the B and R- 
bands than in 2001. While Lll mainly used the 
HST photometry, the flux increase was supported 
by observations of the flux of the [Ca II] A A 7292, 
7324 lines, that will be discussed in this paper. In 
Lll it was proposed that the increase in the op- 
tical flux is caused by reprocessing of the X-ray 
and EUV emission from the ejecta - ring collision. 
This is further strengthened by the changing mor- 
phology of the ejecta as seen in the HST imaging, 
where Larsson et. al. (2012, in prep., in the follow- 
ing L13) find that most of the Ha emission from 
the inner ejecta has a ring-like morphology. This 
is in contrast to the emission in the 1.644 /im [Si 
I]/[Fe II] emi ssion, which main ly comes from the 
core (see also Kiaer et al. 2010() . 

In this paper we discuss ground based high S /N 
observations with high spectral resolution from 
the Very Large Telescope (VLT), complemented 
by STIS observations from HST. A high S/N is 
especially important in order to trace the faint 
line wings to high velocity and also to detect 
additional, w 'eaker lines i i i add ition to the ones 
discussed by ISmith et al.l (|2005l ) and iHeng et al 



(|2006[ ). In particular, we discuss the evolution of 
the reverse shock in time and the evolution of the 
lines from the inner regions of the ejecta. These 
are particularly interesting to monitor since they 
may signal the re-ionization of the inner ejecta by 
the hard radiation from the ring collision. 

In Sect. [2] we discuss the observations and re- 
ductions, and in Sect. [3] we detail the results. In 
Sect. |4]we discuss the implications for the reverse 
shock and for the emission from the ejecta. Sum- 
mary and conclusions are given in Sect. [S] 

2. Observations 

2.1. VLT observations 

The ground based observations were performed 
with the VLT at ESO, using the F0RS2 and 
UVES instruments. SN 1987A has been monitored 
regularly with these instruments since 1999. Table 
ID and Table [2] summarize the UVES and F0RS2 
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observations. Primarily, the purpose has been to 
follow the evolution of the narrow lines from the 
ring collision. We have therefore mainly used the 
UVES high resolution spectrograph. To maximize 
the spectral resolution we have chosen a compar- 
atively narrow slit, 0.8" wide. 

In most seasons, observations with seeing of or- 
der 0.6" were obtained, just enough to spatially 
resolve the ri ng into the northern and southern 
parts (see e.g., Groningsson et al1l2008b ). The po- 
sition angle of the slit was 30° in all observations. 
In Fig. [l]we show the slit orientations and widths 
for the UVES observations at two epochs. 

The F0RS2 observations were obtained with a 
wide 1.6" slit and are accompanied by a local spec- 
trophotometric standard. To estimate slit losses of 
the UVES spectra we have also obtained low reso- 
lution F0RS2 observations with a narrow slit (see 
Table [3]) Additionally, the slit acquisition images 
of the UVES spectrograph have been inspected to 
ensure that the alignment of the instrument slit 
to the supernova was consiste nt throughout the 



However, for the broad Ha and H/3 lines we have 
used the robust 2-D extractions here. 



data set. This, as described in iGroningsson et al 



(|2008bl) . gives us confidence that the error in the 
absolute flux calibration is between not larger than 
30%. Furthermore, the UVES observations em- 
ploy a 12" long slit. This provides a clean obser- 
vation of Ha emission from the LMC at the edges 
of the slit, which it is plausible to assume is in- 
variant over the period of our observations and 
can be expected to be extended with respect to 
the slit (therefore insensitive to seeing). The vari- 
ation of the flux in the LMC lines is well within 
20%. The uncertainties are considered to be al- 
most exclusively due to the position of the super- 
nova with respect to the slit and the seeing rather 
than throughput of the instrument or the atmo- 
sphere. 

The reductions of the F0RS2 data followed 
classical long slit techniques using the FIGARO data 
reduction package. The UVES reductions used the 
ESO pipeline. All wavelength settings were pro- 
cessed using both optimal extraction and a manual 
extraction along the slit. For the Ha and H/3 lines 
the strong emission from the narrow circumstellar 
lines poses a challenge for optimal extraction algo- 
rithms, and the 2-dimensional manual extraction 
and sky subtraction was preferred. The difference 
between the two methods only affects the narrow 
components which do not form part of this work. 




Fig. 1. — Slit position for the UVES observations 
overlaid on HST/WFPC2 F675W images from 
2000 Nov. 13 (left) and 2009 Apr. 29 (right). 
A slit width of 0.8" was used in all observations, 
except for the 1999 observation. 

In this paper we concentrate on the evolution of 
the broad lines from the ejecta and reverse shock. 
To show this more clearly we subtract the narrow 
and intermediate velocity lines by a cubic spline 
interpolation between the extremes of the narrow 
lines. Figure [2] shows the result before and af- 
ter subtraction for Ha for the day 7976 FORS 
spectrum. Even if the resulting line has a smooth 
shape, one should note that especially the peak of 
this line is uncertain because of the narrow line 
contamination. In some cases comparison with 
HST spectra helps in this respect, even if their 
spectral resolution is much lower (Sect. 
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Fig. 2. — Ha profile before (blue) and after (red 
line) subtraction of the narrow and intermediate 
velocity lines for the F0RS2 1200R spectrum (slit 
width 1.62") at 7976 days. 
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Table 1 

VLT/UVES OBSERVATIONS OF SN 1987A 



Epoch 


Date 


Days after 


Setting 


A range 


Slit width 


Resolution 


Exposure 


Seeing 


Air mass 






explosion* 




(nm) 


(arcsec) 


(A/AA) 




(arcsec) 




1 


1999 Oct 16 


4618 


346+580 


303 - 388 


1.0 


40,000 


1,200 


1.0 


1.4 










476 - 684 












2 


2000 Dec 10 - 14 


5040 


346+580 


303 - 388 


0.8 


50,000 


10,200 


0.4-0.8 


1.4 










476 - 684 














Dec 09 - 10 


5038 


390+860 


326 - 445 






9,360 


0.4 


1.4-1.6 










660 - 1060 












3 


2002 Oct 06 - Dec 14 


5738 


346+580 


303 - 388 


0.8 


50,000 


10,200 


0.7-1.0 


1.5-1.6 










476 - 684 














Doc 13 - 14 


5772 


390+564 


326 - 445 






9,360 


0.4-1.1 


1.4-1.5 










458 - 668 














Oct 04 - 05 


5702 


437+860 


373 - 499 






9,360 


0.4-1.1 


1.4-1.5 










660 - 1060 












4 


2005 Mar 21 - Apr 12 


6611 


346+580 


303 - 388 


0.8 


50,000 


9,200 


0.6-0.9 


1.5-1.8 










476 - 684 














Apr 09 - 11 


6621 


437+860 


373 - 499 






4,600 


0.5 


1.6-1.7 










660 - 1060 












5 


2005 Nov 01 


6826 


346+580 


303 - 388 


0.8 


50,000 


2,300 


0.9 


1.4 










476 - 684 














Oct 20 - Nov 15 


6825 


437+860 


373 - 499 






9,200 


0.5-1.0 


1.4-1.6 










660 - 1060 












6 


2006 Oct 01 - 21 


7170 


346+580 


303 - 388 


0.8 


50,000 


9,000 


0.5-0.9 


1.4-1.5 










476 - 684 














Oct 29 - Nov 15 


7196 


437+860 


373 - 499 






9,000 


0.5-1.0 


1.4 










660 - 1060 












7 


2007 Oct 23 - Nov 28 


7565 


346+580 


303 - 388 


0.8 


50,000 


11,250 


0.8-1.4 


1.4-1.6 










476 - 684 














Oct 23 


7547 


437+860 


373 - 499 






9,000 


1.1-1.4 


1.4-1.6 










660 - 1060 












8 


2008 Nov 23 - 2009 Feb 08 


7982 


346+580 


303 - 388 


0.8 


50,000 


9,000 


0.8-1.0 


1.4-1.5 










476 - 684 














2009 Jan 09 - 25 


7998 


437+860 


373 - 499 






11,250 


0.8-1.3 


1.4-1.5 










660 - 1060 












9 


2010 Nov 06 - 15 


8661 


346+580 


303 - 388 


0.8 


50,000 


9,000 


0.7-1.1 


1.4-1.7 










476 - 684 














2010 Oct 19 - Nov 15 


8652 


437+860 


373 - 499 






11,250 


0.6-1.1 


1.4 










660 - 1060 












10 


2011 Nov 06 - Dec 03 


9019 


346+580 


303 - 388 


0.8 


50,000 


9,000 


0.8-1.1 


1.4-1.5 










476 - 684 














2011 Nov 05 - Dec 04 


9035 


437+860 


373 - 499 






9,000 


0.6-1.4 


1.4-1.7 










660 - 1060 













^ Average epoch of spectrum since explosion, 19S7 Feb. 23 . 
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Table 2 

VLT/FORS AND HST/STIS observations of SN 1987A 



Epoch 


Date 




Days after 
explosion 


Setting 


A range 

(nm) 


oiit wiutn 

^circsGC^ 


Resolution 

(A/ L^A) 


Exposure 


Seeing 
(&rcsec) 


Air mass 


STIS 


1999 Feb 21 - 


27 


4384 


750L 


524 - 


636 


0.5 


666 


10, 500 






STIS 


1999 Aug 30 - 


- 31 


4571 


750M 


630 - 


687 


3 X 0.1 


6,000 


3x7, 804 






FORSl 


2002 Dec 30 




5788 


600R 


514 - 


730 


0.70 


1,660 


5,400 


0.7-0.9 


1.4-1.6 


FORSl 


2002 Dec 30 




5788 


600B 


336 - 


576 


0.70 


1,110 


7,200 


0.7-1.0 


1.4 


STIS 


2004 Jul 18 - 


23 


6358 


750L 


524 - 


636 


3 X 0.2 


666 


3x5, 468 






FORS2 


2006 Nov 24 




7213 


600RI 


512 - 


845 


0.7 


620 


1,200 


0.7 


1.5 


FORS2 


2006 Nov 24 




7213 


600RI 


512 - 


845 


1.62 


620 


1,200 


0.7 


1.5 


FORS2 


2006 Dec 21 




7240 


1028Z 


786 - 


962 


1.62 


1,580 


1,320 


0.8 


1.4 


FORS2 


2006 Dec 21 




7240 


1200R 


590 - 


740 


1.62 


1,320 


1,320 


0.9 


1.4 


FORS2 


2006 Dec 21 




7240 


1200R 


590 - 


740 


0.70 


3,060 


1,380 


1.1 


1.4 


FORS2 


2006 Dec 21 




7240 


1400V 


465 - 


596 


1.62 


1,300 


1,380 


0.7 


1.4 


FORS2 


2007 Nov 07 




7561 


600RI 


537 - 


870 


1.62 


620 


600 


0.9-1.1 


1.4 


FORS2 


2007 Nov 07 




7561 


1028Z 


786 - 


962 


1.62 


1,580 


1,740 


1.0-1.2 


1.4-1.5 


FORS2 


2007 Nov 07 




7561 


1200R 


590 - 


740 


1.62 


1,320 


600 


1.3-1.4 


1.4 


FORS2 


2007 Nov 07 




7561 


1200R 


590 - 


740 


0.70 


3,060 


600 


0.9-1.0 


1.4 


FORS2 


2007 Nov 07 




7561 


1400V 


465 - 


596 


1.62 


1,300 


900 


1.0-1.2 


1.4 


FORS2 


2008 Nov 04 - 


- Dec 26 


7950 


600RI 


537 - 


870 


1.62 


620 


1,200 


0.9-1.0 


1.4 


FORS2 


2008 Nov 28 - 


- Dec 03 


7951 


1400V 


465 - 


596 


1.62 


1,300 


1,800 


0.6-1.0 


1.4-1.5 


FORS2 


2008 Dec 26 




7976 


1200R 


590 - 


740 


0.70 


3,060 


600 


0.9 


1.4 


FORS2 


2008 Dec 26 




7976 


1200R 


590 - 


740 


1.62 


1,320 


600 


0.8-1.0 


1.4 


FORS2 


2008 Dec 26 




7976 


1028Z 


786 - 


962 


1.62 


1,580 


1,740 


0.8-0.9 


1.4 


STIS 


2010 Jan 01 




8378 


750L 


524 - 


636 


0.2 


666 


14, 200 






FORS2 


2011 Nov 25 




9040 


600RI 


537 - 


870 


1.62 


620 


1,200 


0.9-1.0 


1.41 


FORS2 


2011 Dec 03 




9049 


1400V 


465 - 


596 


1.62 


1,300 


1,800 


0.6-1.0 


1.55 


FORS2 


2012 Jan 14 




9090 


1200R 


590 - 


740 


1.62 


1,320 


1,400 


0.8-1.0 


1.41 


FORS2 


2012 Jan 14 




9090 


1028Z 


786 - 


962 


1.62 


1,580 


3,480 


0.6-0.9 


1.41 


FORS2 


2012 Jan 23 




9099 


1200R 


590 - 


740 


0.70 


3,060 


1,200 


0.9 


1.55 



^ Average epoch of spectrum since explosion, 1987 Feb. 23. 
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The UVES slit only covers part of the ejecta 
and this affects the flux and the line profiles we 
observe. The major axis of the ring is 1.6" and 
the minor axis 1.1", so only part of the outer 
ejecta, which now fills most of the ring, is within 
the UVES slit (Fig. [J). Because most of the flux 
in the high velocity wings come from the central 
parts of the projected image along the line of sight 
(LOS), these are likely to be less affected than the 
low velocity parts of the line. These come both 
from the supernova core and the outer parts of 
the projected image. Part of the latter emission 
may be missed by using a narrow slit. 

Fig. [3] shows a comparison between one of the 
slit integrated spectra obtained with UVES with 
the 0.8" slit and the F0RS2 spectrum with the 
1.6" slit for Ha. The F0RS2 slit covers most of 
the ejecta, while the UVES slit only covers the cen- 
tral fraction of it. The exact fraction of the total 




-10 -5 5 
Velocity km s"^) 

Fig. 3. — Comparison of Ha taken with F0RS2 
with a 1.6" slit and UVES with a 0.8" slit at 7976 
days and 7982 days, respectively. The lower, blue 
spectrum is the original UVES spectrum, while 
the black spectrum is from F0RS2. The red spec- 
trum is the UVES spectrum multiplied by a factor 
1.86. Most of the remaining difference between the 
UVES and F0RS2 spectra come from velocities 
~ 5000 kms~^, characteristic of the reverse shock 
close to the ring plane and outside the UVES slit. 

flux depends on the seeing and the spatial origin 
of the emission. In Fig. [3] we have subtracted 
the continuum, which mainly comes from the ring 
collision (Sect. [3]). Most of the difference in the 



F0RS2 and UVES line profiles comes from the 
parts of the ejecta outside the slit, i.e., in the NW 
and SE directions. In the ring plane, where most 
of the emission from the reverse shock originates, 
this corresponds to velocities in the range -5000 
- -1-5000 kms~^, while higher velocities, coming 
mainly from material expanding in our direction, 
should fall within the narrow UVES slit. 

In Table [3] we give the ratio of the Ha fluxes 
with the different slit widths at the epochs when 
we have observations with both instruments. This 
indicates that the fraction lost stays nearly con- 
stant with time and is close to proportional to the 
width of the slit. We note that inner ejecta, mov- 
ing across the line of sight at < 4000 kms~^ for 20 
years will still be within 0.33 " of the remnant's 
center. Although slightly spread out by the seeing, 
most of this emission should then fall within the 
0.8" slit of UVES. Most observations were made 
with a seeing < 0.8" and the comparison with the 
FORS observations indicate that the slit losses are 
small for the core component. (See Fig. [H below.) 

2.2. HST STIS Observations 

The HST STIS observations used in this paper 
were carried out in 1999 (day 4381, G750L grat- 
ing, and days 4571-4572, G750M grating), 2004 
(days 6355-6360, G750L grating) and 2010 (day 
8378, G750L grating). The G750L grating has 
a spectral resoltion of ~ 450 kms~^ and covers 
the wavelength interval between 5240 - 10270 A, 
which includes Ha and [Ca II]. The spectral reso- 
lution of the G750M grating is significantly better 
(^ 50 kms^^), but in this observation the wave- 
length coverage was reduced to 6295 - 6867 A, 
which means that only the Ha line is included. 

For the 199 9 observation G750L a wide slit of 
0.5" was used (jMichael et al.ll2003l) at a position 
angle 25.6°. In addition, spectra with the G750M 
grism with three parallel slits of 0.1" width at a 
position a ngle 27° were tak en. For the 2004 ob- 
servation ( Heng et al. 20061) a 0.2" slit with a po- 
sition angle of 0° was placed in three different lo- 
cations, thus covering all of the inner ejecta as 
show n in Fig. [T31 Fina lly, for the 2010 observa- 
tion (jFrance et al.ll2010t ) one 2" slit with posi- 
tion angle 0° was used, covering only the central 
part of the ejecta perpendicular to the slit. Details 
of all observations are summarized in Table [5] and 
in more detail in Table 2 in L13 and we refer the 
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Table 3 

Comparison of Ha fluxes measured with FORS2 and UVES with different slits. 



Days after explosion 
FORS / UVES 


FORS2(1.6")/UVES(0.8") 


UVES(0.8")/FORS(0.7") 


7240 / 7170 


1.9 


0.96 


7561 / 7565 


1.8 


1.18 


7976 / 7982 


1.9 


1.24 


9090 / 9019 


2.0 





reader to that paper for details regarding the data 
reduction. 

3. Results 

In Fig. Hlwe show a subset of the full UVES 
spectra from Dec. 2002 until 2011 Nov., together 
with the STIS spectra from 1999 Feb. and 2004 
Jul. (Sects. 12.21 and l3.2p . Also here we have sub- 
tracted the many narrow lines originating from the 
unshocked and the shocked ring from the UVES 
spectra. These spectra have all bee n de-reddened 
by Er _v = 0.19 mag using the ICardelli et al. 
( 19891 ) reddening law. To show the relative energy 



distribution we plot in this figure A^a . 

Most of the rising continuum towards the UV 
is Paschen and two-photon emission from the ring 
collision, which is apparent when we compare the 
continuum with that in the STIS observations, 
which isolate the ejecta better (Sect. 13. 2p . The 
Balmer continuum below 3646 A from the same 
source is also prominent. 

Superimposed on this are several broad lines, 
in particular Ha, H/3, Mg I] A 4571 and [Ca II] AA 
7292, 7324. We identify an interesting feature at 
~ 9220 A as emission by Mg II AA 9218, 9244. 
In the blue wing of Ha there is most likely also a 
contribution from [O I] AA 6300, 6364. We discuss 
these and some weaker features further in Sect. 

We have fit the total continuum spectrum with 
the sum of free-bound and two-photon continu a 
from H I and He I (e.g.. lErcolano fc Storevll200"6l) . 
The two-photon contribution will be suppressed at 
higher densities, and its relative contribution will 
therefore depend on the density. To include this 
effect we have employed a model atom for H I sim- 



ilar to that used in iKozma fc FranssonI (jl998a[) . 

Our best fits correspond to a temperature of 
~ 3 X 10^ K and a density of ~ 10^ cm~^, although 
the density could be higher without affecting the 
results appreciably. These numbers are reasonable 
for the conditions in the shocked gas of the ring. 

The first thing to note from Fig. |4] is the in- 
crease in the flux of all broad lines during this 
period. A good example is the [Ca II] AA 7292, 
7324 doublet, which shows a steady increase in 
the flux. We also note that the strength of the 
[Ca II] lines in the STIS spectrum are in between 
the UVES 2002 Oct. (day 5702) and 2005 Nov. 
(day 6825) fluxes, which shows that we include 
most of the flux from the central ejecta also in 
the UVES observations, despite the narrow slit. 
There is, however, a large difference in the line 
profile of the H I lines between the UVES and 
STIS observations, with much stronger 'shoulders' 
in the ground based observations. This is caused 
by the larger extraction region along the slit for 
the UVES spectra compared to the STIS spectra 
(see Sect. 13. 2p . This results in a larger contribu- 
tion from the reverse shock in the UVES spectra. 

Because the hydrogen lines and the other lines 
differ substantially in their properties and origin 
we discuss the results of these separately. 

3.1. Hydrogen lines 

Figure [5] shows the evolution of the slit inte- 
grated Ha line profile with time from 2000 Dec. 
to 2011 Nov. from our UVES observations. In ad- 
dition to removing the narrow lines, we have here 
subtracted the continuum level by simply sub- 
tracting a constant flux for each date. Because 
of the slow variation of the continuum level this is 
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Fig. 4. — Compilation of UVES and STIS spectra from the ejecta and reverse shock, with the narrow hnes 
from the ring subtracted. Each spectrmxi has been been shifted upwards by 3 x 10~^^ erg s~^cm~^ relative 
to the previous. The dashed lines give the zero flux level for each date. Both H/3 and H7 are blended with 
Fe I-II lines. The gaps in the UVES spectra are either due to gaps between the grisms or regions severely 
contaminated by lines from the ring. All spectra have been de-reddened with Eb-v ~ 0.19 mag. This figure, 
unlike the others, is displayed in XF\ for clarity. 
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a reasonable approximation, without introducing 
additional parameters. 
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Fig. 5. — Evolution of the Ha profile from 2000 to 
2011 from the 0.8" slit UVES observations. The 
positions of the most important narrow lines from 
the ring collision, which have been blocked out, 
are marked. The continuum, defined by the red 
and blue sides outside of Ha, has been subtracted 
for each spectrum. The blue wing above ^ 8000 
kms~^ is dominated by a blend of [O I] AA 6300, 
6364 from the inner ejecta and the shocked ring. 

Extracting different regions along the slit we 
can separate the northern and southern parts of 
the ejecta as shown in Fig. |6l The northern com- 
ponent is clearly blueshifted, while th e southern is 
redshifted. This was found already bv lSmith et al. 
(j2005l ). and is discussed in detail in L13. We 



have also measured the evolution of the flux in 
the two regions, and find this to be very similar, 
with roughly equal fluxes in the two. 

From Figs. [3] - [5] one can distinguish two 
components t o the line profile, as was done in 



Smith et al.l (|2005[ ). One velocity component 
reaching ~ 4500 kms^^ from the core, and 
one high velocity component reaching velocities 
> 10,000 kms~^. These evolve fairly indepen- 
dently of each other, and to determine the flux 
in each of these components we have assumed the 
broad component to include the part of the line 
profile with velocity higher than ~ 2500 kms^^ 
on the blue and red sides. Between these veloci- 
ties we make a polynomial interpolation between 
the red and blue parts of the profile of the broad 




-10-5 5 10 
Velocity (10^ km s"^) 

Fig. 6. — Comparison of the Ha line profiles for 
the northern and southern parts of Ha at 7982 
days post explosion. 



component, and assume that the flux above this 
represents the contribution from the core and that 
below from the reverse shock (see Fig. [7] for the 
Dec. 2011 (day 9019) spectrum). This is a rea- 
sonable procedure since a radially thin shell, a 
fair approximation to the geometry of the reverse 
shock, should have a flat line proflle (see e.g., 
iMichael et aI1l2003l for simulations with different 
geometries) . 

However, given the limitations of ground based 
observations, there may be substantial systematic 
errors in the fluxes determined. This applies es- 
pecially to the core component, where the inter- 
polation of the line proflle in the red part of the 
line is uncertain. Also the division between the 
reverse shock and the core component introduces 
an additional uncertainty in the core flux. The 
much higher contribution to the flux by the re- 
verse shock implies that it is less affected by this 
procedure. 

In Fig. [8] we show the resulting light curves for 
the core component and the reverse shock compo- 
nent. Note that the UVES sht only covers 0.8" 
of the ejecta and to get the total flux from the 
reverse shock we therefore multiply these by the 
correction factor shown in Table |3l The flux from 
the core component should fall mainly within the 
slit, although part of this may also fall outside 
due to the broadening of the PSF because of the 
seeing. We have in the same way determined the 
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Fig. 7. — Illustration of the flux determination of 
the core component for the Dec. 2011 UVES spec- 
trum. The solid, blue line shows the (smoothed) 
spectrum with the narrow lines removed. The dot- 
ted line in the gaps is the spline interpolation to 
the spectrum in the gaps. The solid, red line is a 
polynomial fit to the reverse shock between —5000 
kms~^ and -1-7000 kms~^, used for the calculation 
of the core flux. The vertical solid lines mark the 
±2500 kms~^ velocity and the dashed, vertical 
line the —1000 kms~^ velocity. The horizontal, 
dashed line shows the continuum level. The flux 
of the reverse shock is calculated between —8000 
kms~^ and -1-10,000 kms^^, with the core com- 
ponent subtracted. 



fluxes from the dates where FORS observations 
are available. This is an important check because 
the 1.6" slit should cover the full ejecta and re- 
verse shock. This is shown as the filled squares in 
Fig. El 

As an additional check of the uncertainty intro- 
duced by the interpolation we also determine the 
core flux in the part of the line between —2500 
kms~^ and —1000 kms^^ not affected by the nar- 
row lines, shown as open circles in Fig. |8] . 

A third way of determining the evolution is to 
plot the monochromatic flux at velocities charac- 
teristic of the central ejecta and reverse shock, re- 
spectively. In Fig. HI we show the flux at -1300 
kms~^, which is as close to the peak as one can 
safely trace Ha without interference with the nar- 
row lines, and at -3000 kms~^, which is at the 
plateau, as well as their difference. We believe 
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Fig. 8. — Upper panel: Evolution of the Ha flux 
from the reverse shock. Filled triangles show the 
flux from UVES, while filled squares are from 
FORS. The reverse shock fluxes from UVES have 
been been multiplied by a factor 2.0 to compen- 
sate for the narrow slit. Lower panel: Same for 
the flux from the core. Filled triangles show the 
UVES flux between ±2500 kms~^, including the 
interpolated part of the line (see Fig. [7|, while 
the open circles are the fluxes from the part of the 
line unaffected by the narrow lines between —2500 
kms~^ and —1000 kms~^. These have been scaled 
by a factor of 4.0 for a comparison with the inter- 
polated total core flux. The red points show the 
corresponding values measured with STIS. 

that, although this only gives the monochromatic 
flux at these velocities, it gives a more accurate 
representation of the flux evolution than the total 
flux. 

For the reverse shock we see that, including the 
correction for the slit width, there is good agree- 
ment between the UVES and FORS fluxes, which 
shows that the more complete UVES set gives a 
good representation of the time evolution. 
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Fig. 9. — Evolution of the monochromatic, contin- 
uum subtracted Ha flux at -1300 kms~^ (squares) 
and at -3000 kms~^ (circles) and their difference 
(triangles) . 



The core fluxes are considerably more uncer- 
tain, both in terms of the absolute flux and time 
evolution. In Sect. 13.21 we estimate the flux from 
the core from high spatial resolution STIS spectra. 
We there find a flux nearly a factor of two higher, 
which gives an estimate of the systematic errors in 
the ground based estimates, mainly caused by the 
contamination by the narrow lines. The fact that 
the total interpolated flux and the directly observ- 
able flux between —2500 to —1000 kms~"'^show a 
very similar increase, however, indicates that the 
relative flux evolution is reasonably accurate. We 
discuss the implications of the light curves further 
in Sect. O 

We can compare the Hct fl ux with that deter- 
mined by ISmith et all (|2005[ ) on 2005 Feb. 25 
(day 6577). On this day they find a total flux 
from the reverse shock component of 1.37(±0.15) x 
10~^^ ergs~^cm~^. This is approximately a fac- 



tor of two higher than what we determine from 
UVES for 2005 Apr. The slit they use is similar 
to ours, 0.8", and the seeing is also similar. The 
position angle is different from ours, P. A. —10°, 
compared to 30°, and may affect the flux some- 
what. However, they calibrate their flux with the 
part of the HST F658N narrowband image of the 
ring that falls within the slit. This should give 
a reasonable calibration, which also compensates 
for the fraction of the light from this region that 



is scattered outside the slit because of the seeing 
in the ground based image. In addition, they add 
a correction from the HST image for the fraction 
of the flux which falls outside of the slit. After the 
correction applied to our data for the UVES slit 
from Table [3l we find a fair ag r eemen t between our 
flux and that of lSmith et al.l ( 20051 ). To compare 
the F0RS2 & UVES data we have applied this 
correction to our light curve. Although this does 
not include all flux, it does result in a consistent 
light curve. The total flux in the reverse shock at 
each time can be found by multiplying the flux in 
Figure [8] by the factor in Tabled 

The maximum velocity of Ha is not well defined 
on either the red or the blue side. On the blue side 
it is difficult to establish the true extent of the 
line due to blending with the broad [O I] AA 6300, 
6364 lines from the ejecta around 9096 kms^^ and 
12,022 kms^-'^, respectively. On the red side nar- 
row and intermediate velocity lines from He I A 
6855 and Fe II A 6872, interfere. We can there- 
fore only establish lower limits to the maximum 
expansion velocities on each side, ^ 9000 kms^^ 
on the blue side while on the red side the line ex- 
tends at least to 11,000 kms~^, where the line is 
still substantially above the 'continuum' level, and 
could extend to considerably higher velocities, up 
to ~ 13,000 kms-i. 

Because of the difficulty in determining the ex- 
act continuum level the evolution of the flux at 
velocities > 10,000 kms~^ is uncertain. At lower 
velocities the flux is increasing monotonously up 
to ~ 7500 days, after which it flattens. The gaps 
in the line profile at high velocities make it diffi- 
cult to establish whether the maximum velocity 
is affected by the interaction with the external 
medium. 

The maximum velocity can be compared to the 
width of the r ed wi ng of Ha at 7.87 years, which 
Chugai et aP (|l997t ) found to extend to - 10, 800 



kms~^. Although the S/N in their observations 
was considerably lower, our velocity measur e ments 
are consistent with theirs. iMichael et al. (l2003h 
report diffuse Ha emission out to ^ ±12,000 
kms^^ from the central direction of the ejecta 
in their Oct. 1999 observation, presumably the 
same emission as is giving rise to the high velocity 
wing in our observations. This is also supported 
by the fact that this part of the line is similar in 
the northern and southern part of the spectrum 
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(Fig. 16]), as expected from a direction in the pro- 
jected center. As we show in Sect. 13.21 also the 
STIS 2004 and 2010 spectra give a velocity close 
to this. 

Most of the change in flux of Ha takes place 
at low and intermediate velocities. A clear break 
in the profile can be seen at ±4000 km s^^ (Fig. 
[5|). Above this velocity the line profile shows a 
gradual decrease to the maximum velocity. The 
main change in flux of the UVES profiles occurs 
below ~ 6000 kms~^ for both the blue and red 
wings, as can be seen in Figs. [5] and [9] At higher 
velocities the profile is nearly unchanged. 

This can be investigated in somewhat more de- 
tail from the FORS spectra in Fig. (TU] When 
comparing the spectra from days 7240 and 7976 
we see that the increase in flux only occurs for ve- 
locities below ^ 5000 kms~^. If we now compare 
the spectra from days 7976 and 9090 the change 
only occurs for velocities below ~ 4600 kms~^. 
Although these velocities are uncertain, and de- 
pend on e.g., the exact continuum level, we there- 
fore note that the velocity below which the flux 
increases is getting smaller with time. 

Below ~ 2500 kms~^ there is a linear rise of 
the profile with time, which is likely to be caused 
by a superposition of the core component and the 
flat reverse shock profile. This can clearly be seen 
in the STIS observations discussed in Sect. 13.21 
On the red side no flat section is seen, which 
is explained by the fact that the region between 
3300 — 6000 kms~^ is blocked by lines from the 
ring. We also note that there is an asymmetry be- 
tween the red and blue sides, with a considerably 
lower flux on the red (southern) side. 

When comparing the line profiles from the 
southern and northern parts (Fig. ^ it is clear 
that the blue wing of the 'flat' part of the line orig- 
inates from the northern half of the ejecta/ring, 
while the opposite is the case for the red part, as 
is expected from the geometry. The higher flux 
from the NE side of the debris, which produces the 
blue-shifted side of the line profile, is consistent 
with the stron ger e.iecta-ring inter action seen on 
that side fe.g.. iRacusin et al. 



20091) . 



Figure [TT] shows H/3 together with Ha for days 
7982 - 7998. While the continuum level of Ha is 
well determined, H/3 is in a crowded region and the 
continuum is more difficult to define, which intro- 
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Fig. 10. — Comparison of Ha profiles from F0RS2 
with a 1.6 " slit from 2006 Dec. (day 7240) to 
2012 Jan. (day 9090) (bottom to top). The upper 
panel shows the original spectra, while the lower 
plot shows spectra with the narrow lines removed 
and the continuum subtracted. Note that it is 
mainly the part of the line core below ^ 5000 
kms~^ which has increased in flux during this pe- 
riod. The position of the [O I] AA 6300, 6364 dou- 
blet from the core, which influences the blue wing 
of the Ha line, has been marked. 



duces some uncertainty. By scaling the H/3 flux 
by a factor 3.0, as in Fig. [TTl we find good agree- 
ment between the two lines for the blue wing. The 
red wing of H/3 is, however, considerably stronger 
compared to Ha. This is most likely a result of 
blen ding with othe r lines. In their 1995 observa- 
tions IChugai et al.l (|l997l) identify two Fe II lines 
at 4889 A and 5018 A. In addition, the 'bump' at 
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^ 4000 kms ^ can be identified with the strong 
Fe II a 6S - z 6P° A 4924 transition. The ve- 
locities of these hnes are marked in Fig. [TTJ We 
also note the lower velocity of the blue wing of II/3 
compared to Ha, which confirms that this part of 
'Ha' is caused by blending with the [O I] AA 6300, 
6364 lines. The lines therefore complement each 
other in filling several of the 'gaps' caused by the 
lines from the ring collision, as well as blending 
with other ejecta lines. 




Velocity (10 km s 



Fig. 11. — Comparison of the Ha (blue) and H/3 
(red) UVES profiles for day 7982 - 7998. The 
flux has been scaled by a factor 3.0 compared to 
Ha. The red wing is probably blended with three 
Fe II lines, marked with their positions in velocity. 



Including a correction for reddening, the ob- 
served Ha/H/? ratio of 3.0 corresponds to Ha/H/3 = 
2.5. There is no indication that this ratio varies 
over the line profile. The measured ratio is, how- 
ever, sensitive to the chosen level of the con- 
tinuum, which is especially problematic for H/3, 
where there are few regions free from other lines 
(see Fig. 1111) . Given this unc ertainty, our Une ra - 
tio is similar to that found bv lChueai et al.l(ll997h 
at 7.87 years, Ha/H/3 = 3.8. At that epoch the 
flux was, however, dominated by the inner ejecta, 
with little contribution from the reverse shock. 
Because the emission from the reverse shock and 
from the ejecta are produced by different mech- 
anisms, there is no reason to expect these ratios 
to remain the same (Sect. 14. ip . This ratio can be 
compared to modeling of t he Balmer emis s ion in 
supernova remnants, where [Chevalier et al. ( 1980[ ) 



find an Ha/H/? ratio of 3 - 5 depending on the 
optical depth in the Lyman lines. The uncertainty 
in the value determined here is large enough that 
it is unclear if there is any diagnostic value in this 
ratio. 

In addition to these lines one can in Fig. |4] 
identify a strong line coinciding with H7, which 
is clearly blend ed with other l ines. Based on the 
obser vations in IChuga.i et al.l (|l997() and model- 



mg m lJerkstrand et al.l (|201lr ) these are mainly Fe 
I, and in some cases als o Fe II, lines. The Ca 
I A 4226 line, found by Ijerkstrand et al.l (|201ll ) 

to be strong in the 8 year spectrum, may also be 
present. Higher members of the Balmer series may 
contribute below ^ 4200 A. 

3.2. Comparison with STIS observations 
of Ha 

HST narrow slit spectroscopy of Lya and Ha 
between 1 999 - 2004 ha s been discu ssed exten- 
sively by iMichael et al. ( 2003 ) and Heng et al. 
( 20061 ). and recent observations from 2010 have 



been discussed by lFrance et al. ( 2010l ). These pa- 
pers were focused on the reverse shock properties, 
while we here emphasize the core component and 
the connection to our ground-based observations. 
Although the spectral resolution of the G750L 
grating is only ^ 450 kms~^ (FWHM), which is 
more than one order of magnitude less than UVES, 
it is adequate for the broad lines. The most im- 
portant advantage of the STIS data is that it is 
possible to directly relate different spatial regions 
to their velocities. In the following discussion we 
will focus on the observations from 1999 and 2004, 
where the slits covered the full central region of 
the ejecta. The 2010 observation only covered the 
central 0.2" and is therefore difficult to compare 
with the other observations. 

The last ejecta spectra not strongly affected by 
the ring collision are the STIS G750L and G750M 
spectra from 1999. The full G750L spectrum from 
1999 Feb. (day 4381) is shown in Fig. gl This 
spectrum was taken with a 0.5" slit, which covered 
the full inner ejecta at this epoch. The spectra 
taken with the G750M grating (days 4571-4572) 
only cover the Ha line, but offer a much better 
resolution of 50 kms^^. The observations were 
performed with three 0.1" slits, which together 
cover most of the inner ejecta. In the upper panel 
of Fig. [12] we show a comparison of the Ha profiles 
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from the G750L and G750M gratings. Both spec- 
tra were extracted from the central ±0.3" along 
the slits. There is clearly a good agreement be- 
tween the two spectra, with the only notable dif- 
ference being that the medium resolution profile 
contains narrow lines from Ha, [N II] AA 6548, 
6584 due to contamination from the outer ring, 
passing through the ejecta. We will discuss the 
line profile more below. 

In the later 2004 spectrum, the ring interaction 
is much more prominent. Figure [T3l shows an ACS 
image of SN 1987A taken in the F625W fiher on 
2003 Nov. 28 (day 6122) with the three slit posi- 
tions, each 0.2" wide, as well as the 2-dimensional 
spectrum from t he central sli t . The latter is sim- 
ilar to Fig. 4 in iHeng et all (|2006l) . We include 



it here to illustrate the area of extraction for the 
spectra. 

The full STIS spectrum from the central area 
indicated in Fig.[T3]is shown in Fig. 01 The degree 
of contamination from the outer ring and scat- 
tered light in the central part of the line is difficult 
to estimate. We plot the 2004 spectrum together 
with the 1999 low resolution spectrum in the lower 
panel of Fig. [T2] The 2004 spectrum is the to- 
tal over all three slits, scaled to the same level as 
the 1999 spectrum outside of ±2000 kms"^ We 
first note the similar shape of the wings between 
2000 — 5000 kms~^ on both the red and blue sides. 
There is little evolution between these epochs in 
these velocity intervals. The central part of the 
2004 spectrum is dominated by scattered light 
from the shocked ring emission as evidenced by 
the fact that the Ha / [N II] ratio is similar to t hat 
of the shocked ring ( Groningsson et aDl2008a[ ). 

The 2-dimensional STIS spectra clearly display 
a high velocity component close to the ring on both 
the red and blue sides. This is the freely expand- 
ing, neutral reverse shock component. Most of this 
emission comes from ~ 3000 — 8000 kms~^ on the 
blue side and from ~ 4000 — 7000 kms~^ on the 
red side. The blue component is fairly uniformly 
distributed over the three slit positions close to 
the northern part of the ring, while the fainter 
component on the red side mai nly comes fr o m the 

(l2006h 



Heng et al 



south-eastern side of the ring, 
find that the flux from the reverse shock compo- 
nent in Ha increased by a factor 2 — 3 from 1999 
Sep. 18 to 2004 Jul. 18 (days 4589 and 6355, 
respectively). 




Fig. 12.— Upper panel: The STIS G750L (blue) 
and G750M (red) Ha line profiles from 1999 Feb. 
(day 4381) and Aug. 1999 (day 4571-4572). The 
G750M spectrum has been multiplied by a factor 
1.5 to compensate for the smaller total slit width 
(3 X 0.1 ") compared to the G750L spectrum (0.5 
"). Lower panel: STIS G750L Ha line profiles 
from 1999 Feb. (day 4381) and 2004 Jul. (days 
6355-6360). The 2004 spectrum has here been 
multiplied by a factor 1.6 for comparison with the 
1999 spectrum. The region between ±1500 kms~^ 
in the 2004 spectrum is dominated by scattered 
Ha, [N II] AA 6548.1, 6583.5 emission from the 
ring collision. 



The increase in the northern component seen by 
Heng et al. ( 20061 ) can also be seen in our UVES 
and FORS spectra in Fig. [5] on the blue side at 
^ —5000 kms^^. A similar increase on the red 
side from the southern part of the reverse shock 
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Fig. 13.— Left panel: The HST/ACS F625W image from 2003 Nov. 28 (day 6122). The vertical black lines 
show the three slits positions used for the 2004 Jul. (days 6355-6360) STIS observations. The slits are 0.2" 
wide. North is up and west is right. Right panel: The spectrum of the central slit from the STIS observations, 
taken with the G750L grating. The spectrum is centered on zero velocity for Ha. The horizontal red and 
black lines in both panels show the regions that were used to produce the spectra in Fig. [141 The spectrum 
in the right panel has been stretched in the y-direction to match the spatial scale of the image on the left. 
The velocity scale is in 1000 kms~^ from the rest wavelength of Ha. 
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Fig. 14.— STIS G750L spectra from 2004 Jul. 18- 
23 (day 6355-6360). The spectra were extracted 
from the areas shown in Fig. [131 and are shown 
with the same color coding as in that figure. The 
red, bottom spectrum was extracted from the cen- 
tral ±0.3" and the black, top spectrum is from 
±0.9". 

can also be seen. The high spatial resolution ob- 
servations by HST are therefore consistent with 
our VLT observations. 

We extracted one-dimensional spectra from the 
regions marked by red and black lines in Fig. 1131 
and added the resulting spectra from all three slit 
positions. Figure [14] shows the resulting Ha line 
profiles with different sampling in the N - S direc- 
tion from the center. The similarity of the Ha line 
from the full area (area within the black lines in 
Fig. [131) with the UVES 2005 Mar. spectrum in 
Fig. \E\is noteworthy. Although the S/N, as well 
as the spectral resolution, of the STIS spectrum is 
much worse than the UVES spectrum, one can rec- 
ognize several similarities. In particular, the fiat 
regions at 3000 - 8000 kms~^ on both the red and 
blue sides are apparent in this integration. The 
low velocity part of the line is completely domi- 
nated by the narrow lines from the ring collision. 

In the next extraction from the inner region be- 
tween -0.3" and +0.3" (red lines) the high veloc- 
ity component of the line decreases by a factor 
> 3. This is clear evidence that this spectral fea- 
ture originates in the reverse shock, appearing as 
a thin streak in Fig. [T31 (right panel). The cen- 
tral region is dominated by the ejecta component. 



There is, however, still a flat region at both high 
positive and negative velocities, originating from 
the reverse shock in the line of sight towards the 
center of the ejecta. 

The spatial resolution of STIS allows us to iso- 
late the core component, and in this way check 
the VLT observations. Assuming that the emis- 
sion from the core has the same profile as in 1999 
(Fig. [T^ we find for the Ha flux on day 6355 
that the flux is nearly a factor of two larger than 
the fiux from the UVES measurement that is clos- 
est in time (from day 6622) (Fig. [H). A com- 
parison between Fig. [5l and the red spectrum in 
Fig. [m shows that there are two main reasons 
for this. Firstly, the blue wing of the UVES core 
component is lost due to the interpolation between 
— 1000 and 2000 kms^^and, secondly, some of the 
UVES low-velocity component (especially on the 
red side) has been lost in the process of removing 
the narrow lines. 

The line profile from the low resolution STIS 
spectra of the ejecta shows some indication of an 
asymmetric velocity line profile, with more emis- 
sion on the blue side. Unfortunately, the low ve- 
locity part of the line profiles in the G750L spec- 
tra is severely affected by scattered ring emission. 
We therefore concentrate on the medium resolu- 
tion spectrum from 1999. To study the asymme- 
try of the line we have reflected the line profile 
around zero velocity, and compared this with the 
original in Fig. [151 

When we compare the red and blue line profiles 
we see that for > 2000km s~^ the blue and red 
wings are nearly identical. At lower velocities the 
blue wing is ^ 15% stronger. We therefore find 
that solely on the integrated line profile there is 
only weak evidence for an asymmetry from the 
Ha line. 

The 2D spectra in Fig. [131 on the other hand, 
show strong blue emission from especially the 
northern side, while there is little emission from 
the southern side. This originates in the inner 
ejecta and extends from zero velocity to ~ —5000 
kms~^. This red/blue asymmetry is consistent 
with dust absorption of the emission from the far 
side of the ejecta. 

The message we get from the integrated line 
profiles and the 2D spectra may seem somewhat 
contradictory. It does, however, show that even 
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Fig. 15. — Comparison of the red and blue wings 
of Ha from the STIS G750M 1999 spectrum. The 
blue spectrum is the original spectrum, while the 
red has been reflected around zero velocity. The 
narrow Ha and [N H] A A 6548, 6584 lines seen in 
Fig. [T^have been removed . 



strong asymmetries may cancel in the integrated 
spectra and ideally one needs the full 2D (or even 
better 3D) information. This is discussed in more 
detail in L13. 

3.3. Metal lines from the inner core 

Although not covering the blue part, the STIS 
1999 (day 4384) spectrum in Fig. |4]shows the most 
important metal lines, as well as Ha. The distri- 
bution of the different elements can be inferred 
from their line profiles. In Fig. [16] we show these, 
centered on the blue doublet components. When 
we compare the blue wings of these lines we see 
that the Na I and [Ca II] lines have very similar 
profiles to Ha, while the [O I] A 6300 line has a 
smaller b lue extent. This is con sistent with the 
results in Ijerkstrand et al.l (|201l[ ) where the Na I 
and [Ca II] lines are dominated by emission in the 
H envelope, while the [O I] lines arise mainly in 
the O rich core. 

Note that the [O I] A 6364 line sits on top of 
the flat reverse shock component from the Ha line 
(Fig. [T^ . which partly explains its comparable 
strength to the [O I] A 6300 component (usually 
1/3 of the latter). In addition, there may be a 
contribution to these lines from Fe I A A 6280, 6350 
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(see Fig. 4 in lJerkstrand et al.ll201lf ). 



Fig. 16. — Comparison of the Ha, Na I AA 5890, 
5896, [O I] AA 6300, 6364, and [Ca II] AA 7292, 
7324 from the STIS G750L observation from 1999 
Feb. (day 4384). For a direct comparison of the 
different profiles the fluxes of the lines have been 
scaled by a factor given in the figure. The doublet 
lines have been centered on the blue components. 



Except for Ha, easily identified broad lines in 
the UVES spectra are [Ca II] AA 7292, 7324 and 
Mg I] A 4571. There is also a clear feature coin- 
ciding with the [O I] AA 6300, 6364 doublet. In 
addition. Fig. |4] shows an increasingly strong line 
feature at ^ 9210 A. In Sect. 13.41 we identify this 
with Mg II AA 9218, 9244. Of these the [Ca II] 
lines, and at the last epochs the Mg II lines, are 
both the strongest and best defined. 

As mentioned in the introduction, and dis- 
cussed in Section 14. 3[ the energy input after ~ 
5000 days comes from X-rays produced by the col- 
lision with the ring more than from radioactivity. 
To compare the line profiles at the latest stages 
when X-rays dominate, we show in Fig. [17] Ha, 
[Ca II] AA 7292, 7324 and the Mg II AA 9218, 9244 
lines from day 9019 - 9935. We have here removed 
the narrow lines, as before, and also subtracted 
a continuum level from each line. To more eas- 
ily compare the core component of the lines we 
show in the lower panel the same lines normal- 
ized to the same level for the blue wing, which 
is well defined for all three lines. To isolate the 
core component we have for Ha put the 'contin- 
uum' level close to the 'shoulder' on the red wing 
(at ~ 3.5 X 10~^^ ergs~^ cm~^ in the upper 
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panel). 

When comparing the [Ca II] and Ha Hnes in 
Fig. [17] the former lacks the flat, boxy part of the 
line profile. As the lower panel shows, the central 
part of Ha and the blue wing of the two other 
lines, below 4000 — 5000 kms~^ are, however, 
very similar and argue for a similar origin. The 
full extent on the blue side is difficult to estimate 
because of contaminating lines from the ring, but 
reaches at least 4000 kms^^. On the red side the 
lines are less well-defined, and also affected by the 
blue components of the [Ca II] and Mg II doublets. 
The Mg I] A 4571 line has a lower S/N, but its 
profile is consistent with that of the [Ca II] line. 
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Fig. 17. — Upper panel: Comparison of the Ha, 
[Ca II] and Mg II AA 9218, 9244 profiles at 9019 - 
9035 days. Lower panel: Normalized line profiles 
of the same lines. For Ha the 'continuum' level 
has been set at ^ 3.5 x 10^^^ ergs^^ cm^^ A^^ in 
the upper panel. Note the similar profiles of the 
blue wings of the [Ca II] and Mg II lines to the 
core of Ha. 

Similarly to the other broad lines the line pro- 
file of the [Ca II] line is contaminated by several 




-10 -5 5 10 



Velocity (itf km s"') 

Fig. 18.— Same as Fig. [5] for the [Ca II] AA 
7292, 7324 lines. Velocities are measured from 
the A 7292 component. For comparison we have 
also added the STIS spectrum for day 6358 as the 
dashed line. 

strong lines from the ring collision. We can never- 
theless trace most of it by interpolation (Fig. [THl) . 
This provides a good representation of the true 
line profile as confirmed by the comparison of the 
[Ca II] line profiles from the Jul. 2004 STIS spec- 
trum and the 2004 Apr. UVES spectrum, shown 
as the dashed line in the figure. 

Figure [18] shows that the general line profile re- 
mains roughly the same over the whole period. In 
particular, the maximum velocities of the red and 
blue wings remain at ~ ±5000 kms^^. This and 
the very different shape of the line in comparison 
with Ha argue strongly for this line originating 
in the core and not from the reverse shock re- 
gion. This is also confirmed by the spatial location 
in the STIS spectrum of the [Ca II] region (Fig. 
I19p which shows a central component extending to 
~ 5000 kms-i. 

The absence of a reverse shock component in 
the [Ca II] lines (as well as other ionized species) 
is mainly a result of the low abundance in the 
envelope. The combination of a forbidden tran- 
sition and a low ionization potential also means 
that the number of excitations per ionization will 
be low. This is in contrast to the Li-like ions, 
like C IV, N V and O VI, which have resonance 
lines far below the ionization level. This results 
in a large number of excitations for each ioniza- 
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tion ( Laming et al. 1996f) . This compensates for 
the low abundance and produces hues that are as 
strong as those of H and He. This is, however, not 
the case for the [Ca II] and the other forbidden 
hues in the optical, like the [O I] and Mg I] lines. 
We also note that because of the positive charge 
the Ca II ions will, in contrast to H I, be affected 
by the magnetic field in the shock. Their velocity 
will therefore quickly be isotropised and the ve- 
locity broadening will therefore be characteristic 
of the post-shock velocity. 
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Fig. 19. — Two dimensional STIS spectra from 
2004 (lower) and 2010 (upper) of the Ha (left 
strong line) and [Ca II] (right fainter line) regions 
from the central slit in Fig. [131 To show the faint 
[Ca II] line better the intensity levels have been 
shifted compared to Fig. [T31 Note the absence of 
any high velocity reverse shock component similar 
to that of Ha for [Ca II] . We also show the posi- 
tion of the [O I] doublet, which can be seen also 
in the ejecta with a velocity distribution similar to 
the [Ca 11] lines. The velocity scales are in 1000 
kms^^ from Ha and [Ca II], respectively. 

The upper panel of Fig. [20] shows the evolu- 
tion of the total fiux of the [Ca II] lines from the 
UVES observations. The squares show the flux de- 
termined from the interpolated line profile where 
the missing sections of the line in Fig. [TH| have 
been replaced by spline interpolations. To esti- 
mate the uncertainty introduced by this interpo- 
lation we have also calculated the fiux from the 
line, omitting these regions. These are shown as 
triangles. For a comparison with the interpolated 
line we have multiplied these fiuxes by a constant 
factor 2.5. As an extra check we also show the fiux 



from the STIS 2004 observation measured from 
the central area shown in Fig. [TSj This fiux 
(7.4 X 10~^^ erg cm~^ s~^ on day 6355 after the 
explosion) is similar to that found from UVES, 
and well within the error bars. 

We see in Fig. [20] that the fiux increased mono- 
tonically by a factor of 4 — 6 up to day ^ 7000. 
At later epochs it is nearly constant. This increase 
in fiux is an important indication that there is ad- 
ditional energy input to the inner ejecta in addi- 
tion to the radioactive energy source. In Lll these 
lines were used as an independent confirmation of 
the HST photometry, where the spatial informa- 
tion was combined with the photometry to show 
that the increased flux was coming from the inner 
ejecta, and not from the reverse shock region. 

At these stages, the [Ca II] AA 7292, 7324 
lines are mainly excited by fluorescenc e of UV 
emission through the H fc K lines (iLi et al 



1993; Kozma fc Franssonll 1998b : Jerkstrand et al 



20111) . Therefore the Ca II triplet AA 8498, 8542, 
8662 is expected to have the same total fiux as the 
AA 7292, 7324 lines, unless collisional de-excitation 
is important. The STIS spectrum in Fig. [4| shows 
that there indeed is a line feature at this wave- 
length range, and with a fiux consistent with that 
of the AA 7292, 7324 lines, although it is more 
smeared out in wavelength than the A A 7292, 7324 
lines. 

We also note the presence of an additional 
broad line with a peak wavelength close to that 



of He I A 5876/Na I AA 5890, 5896. IChueai et al 
(119971 ) claim that this is Na I and based on mod- 
eling by [jerkstrand et al. ( 2011 ) this is a likely 
identification. 

3.4. The A 9220 A feature 

The broad line at ~ 9220 A has no obvious 
identiflcation. The full extent is ±3500 kms~^, 
indicating an origin in the core of the supernova. 
It is unlikely that this is a blend of lines from the 
shocked ring, which have a width of ±300 kms""'^, 
since the individual lines would then be resolved. 
This is also confirmed from a direct inspection of 
the STIS spectra, which shows that the A 9220 
feature has the same spatial distribution as the 
[Call] AA 7292, 7324 hues. Possible identifications 
include the Paschen 9 line at 9229.0 A. In this 
case we would, however, expect to see PIO at 9015 
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Fig. 20. — Upper panel: Evolution of the [Ca 
II] A A 7292, 7324 flux with time with UVES and 
STIS. The squares are from UVES and the circle 
point at 6355 days from STIS. The triangles show 
the flux from the regions of the line not affected 
by the narrow lines, multiplied by a factor 2.5. 
Lower panel: Same for the Mg II AA 9218, 9244. 
The fluxes from the unblocked regions of the lines 
have in this case ben scaled by a factor 1.4. 

A nearly as strong and P8 9546 a factor ^ 1.5 
stronger. Neither of these are seen at this level. 
In addition, for pure Case B recombination this 
line is expected to have a flux of ~ 3% of H/3 
(jHummer h, Storevi Il987l) , which is too weak to 
explain the line. 

Another possible identification is a blend of Fe 
II lines, powered by Lya fluorescence. This is the 
mechanism which produces the Fe II lines in SN 
1995N , but in that case, the Ly g was produced 
in the circumstellar interaction (Fransson et al 



20021 ). 



Fe II spectra including Lya excitation from the 
first excited state, a*D, to higher levels, in partic- 
ular u'^P, u'^D and v'^F. The cascade from these 
to lower levels result in a number of strong Fe II 
lines, with the strongest being the 9122.9, 9132.4, 
9175.9, 9178.1, 9196 .9, 9203.1 and 9 204.6. 

Observationally, iHamann et al.l (|l994l ) see a 



strong feature at 9216.5 in the spectrum of Eta 
Carinae. They identify this mainly with Fe II AA 
9203, 9204. Based on their observation, also Fe II 
8490.1 is then expected to be strong. These au- 
thors claim all these lines to be pumped by Lya. 
A weaker broad feature at this position is prob- 
ably present in the STIS spectrum in Fig. HI so 
this may be consistent. The same fe ature is also 
seen i n the symbiotic nova PU Vulbv fRudv et al 



Sigut fc PradhanI (|l998l . |2003[) have calculated 



()l999l) and again identified as Fe II A A 9176-9205 
For the Lya pumping mechanism to work the 
Lya flux needs to be strong and the population of 
the Fe II a^'D level large enough for fluorescence to 
be efficient. In addition, the width of the Lya line 
should be large enough to have a substantial flux 
at the Fe II line. The first and third requirements 
are certainly fulfilled. Both the Lya flux from 
the ring collision fe.g.. llleng et al.ll2006i and refer- 
ences therein ) and from the radioac tive excitation 
of the ejecta ( Jerkstrand et al.lboil ) are strong. A 
major problem is, however, to understand the exci- 
tation to the a^D level from which the fluorescence 
takes place. The excitation energy corresponds to 
11, 000 K, which is much larger than the temper- 
ature expect ed in the radioactively powered core, 
100 - 200 K (jJerkstrand et al.ll201l[ ). The level is 
powered by non-thermal excitations and recombi- 
nations, but these are unlikely to be sufficient to 
make the line optically thick. We therefore con- 
clude that the Lya/ Fe II fiuorescence is unlikely 
to work here. 

A more plausible identification of the A 9220 
feature is Mg II AA 9218, 9244 from the 4p 2p° 
level. This line can arise either as a result of 
fiuorescence by Lya or Ly/3. The 5p ^P° level 
is connected to the Mg II ground state by AA 
1026.0, 1026.1, in nearly perfect resonance with 
Ly/3 at 1025.7 A. This decays to 4p ^P" via either 
the 5s ^S or 4d level, with emission at 8214, 
8235 and 7877, 7896, respectively. This has ear- 
lier been suggested to ex plain UV and IR lines seen 
in QSOs and Seyfert s bvlGrandi fc Phillii^ (Il978h 
and lMorris fc WardI (|l989l ). In our case this has a 
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problem in that the AA 7877, 7896 and AA 8214, 
8234 Hnes are expected to have similar strengths 
as the AA 9218, 9244 hnes, but are not seen in the 
spectra. The efficiency of the Ly/3 / Mg II flu- 
orescence compared to branching into Ha is also 
expected to be low. 

Instead, we propose that the lines arise as a 
result of Lya fluorescence directly to the 4p 
level. The wavelength differences to the multi- 
plet levels are 1239.9, 1240.4 A. A velocity shift of 
~ 6000 kms^^ will therefore redshift the Lya pho- 
tons into these lines, unless other optically thick 
lines are present at shorter wavelength. A veloc- 
ity difference of this magnitude can e.g., arise from 
a photon emitted on the far side of the core and 
absorbed on the near side, 

In the analysis of the 8 year spectrum in 
Jerkstrand et al.l (j201lh it was indeed found that 
the Mg II line dominates the production of the A 
9220 emission. Magnesium is mainly in the form 
of Mg II in the hydrogen rich regions, and the op- 
tical depth of the AA 1239.9, 1240.4 lines are larger 
than unity in the central regions within the super- 
nova core, i.e. within 3000 — 4000 kms^^. This 
remains true also at ^ 20 years and we therefore 
expect the fluorescence to be effective also at these 
epochs. We also find that the lines at AA 10,914 
-10,952 should have a strength a factor ~ 2 fainter 
than the 9220 feature. Inspecting on e of the few 



existin g NIR spectra of the ejecta in iKiaer et al 
( 2010l ) there is a strong feature at this wavelength, 
although this region is noisy and contaminated by 
emission from the ring. This therefore supports 
the identification above. Other lines expected are 
the AA 2929 , 2937 hnes in the UV. These are, how- 
ever, likely to be scattered by the many optically 
thick resonance lines in the UV. 

As is apparent from Fig. [H the flux of this line 
increases rapidly with time. In Fig. [501 we show 
the total line flux from the spline interpolated line 
profile as squares. As in the previous section, we 
also show the flux from the unblocked sections of 
the line as triangles. Because a smaller fraction of 
the line is affected by narrow lines, we only have 
to multiply this by a factor of 1.4 to bring it to 
the same level as the full, interpolated flux. The 
first point at 5038 days shows a negative flux due 
to a combination of the noise and level of the sub- 
tracted continuum. It is therefore an estimate of 
the errors in the flux determination. 



When we compare it with the evolution of the 
[Ca II] lines in the upper panel we see the same 
basic evolution. The scatter of the points is, how- 
ever, smaller due to the smaller fraction of line 
blocked by narrow lines. 

4. Discussion 

The broad lines extracted from the spectra of 
SN 1987A in this paper represent observations 
of a supernova ejecta in the transition from ra- 
dioactively powered supernova to one powered by 
the interaction with the CS environment. Only a 
few supernovae dominated by circumstellar inter- 
action already after the first year, like SN 1979C 
(|Milisavlievic et al. I l2009l) . can compete with this. 
Compared to these we have in the case of SN 
1987A also spatial information, which is crucial 
in interpreting the observations. 

The emission from the ejecta can clearly be 
separated into two components. One of these 
comes from the inner regions of the supernova 
ejecta, well represented by the [Ca II] AA 7292, 
7324 lines. This is the same component as was 
see n in the spec t rum a t 5 - 6 years and 7 . 8 year s 
by IWang et all (|l996l ) and IChugai et all (| 19971 ). 
respectively. The other component, seen in the 
Balmer lines, comes from the reverse shock, re- 
sulting from the interaction with the circumstel- 
lar ring. The Balmer lines have, however, also a 
lower velocity component from the inner parts of 
the ejecta. Because of the different origins of the 
two components we discuss them one by one in the 
next sections. 



4.1. Reverse shock 

The reverse shock 
discussed b y e.g.. 



in SN 1987A has been 



Michael et al. 


(1998b. 20031: 


Heng et al. ( 


20061) and France et al. 



(l201fl[ ). The most important results in this paper 
regarding the reverse shock are the evolution of 
the line profiles, as well as fiuxes. 

As mentioned in Sect. [1] ISmith et al 



(|2005l ) 



discuss the evolution of the fiux from the reverse 
shock extensively. A particularly interesting point 
was their prediction that, around 2012-14, the Ha 
fiux from the reverse shock should reach a maxi- 
mum and then gradually vanish. The rationale is 
that the EUV and soft X-ray flux from the shocks 
ionizes the neutral hydrogen before it reaches the 
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reverse shock. The protons are then accelerated 
in the shock and do not emit any hne radiation, 
except by charge transfer with any unshocked H I 
still present. 

From the Ha luminosity and the estimate that 
0.2 Ha p hotons are emitted for each hydrogen 



ionization Smith et al 



estimate a flux of 
8.9 X 10^^ H I atoms crossing the shock per 
second at an age of 18 years, corresponding to a 
density of n(HI) -- 60 - 70 cm"^. From the X- 
ray luminosity they estimate an ionizing luminos- 
ity of 3.7 X 10^^ s~^ by number, corresponding to 
^ 8 X 10"^^ ergs~^. They therefore estimate that 
~ 4% of the incoming hydrogen atoms become ion- 
ized as they cross the shock. 

The observations by ISmith et al. ( 2005[ ) only 
extended to 2005 Feb. (day 6577), and had rather 
large observational uncertainties arising from the 
different instruments and lines used, limiting the 
extent to which firm conclusi ons with regard to 
this prediction could be made. Heng et al. ( 2006[ ) 
used available STIS observations and showed that 
Ha had increased by a factor of ~ 4 from 1997 
to 2004. Unfortunately, there was little spatial 
overlap between the various epochs and the un- 
certainties in the flux were therefore large. 

In spite of not covering the full ejecta because 
of the narrow slit, our UVES observations in Fig. 
[5] have the advantages of being obtained with the 
same instrumental setup, as well as monitoring the 
same area of the ejecta and ring. In addition, our 
time coverage is a factor of two longer. Further, 
the F0RS2 observations in Fig. [10] cover the full 
ring (except for minor slit losses due to the seeing) . 
These provide a good estimate of the total flux, as 
well as a check of the evolution as determined from 
UVES. 

The light curve of the reverse shock (Fig. |8]) 
exhibits a steady increase in the flux by a fac- 
tor of ~ 3.5 between days 5000 and day 7500. 
However, after ^--7500 days there is a clear flat- 
tening of the light curve. The total increase from 
5000—9000 days is a factor 4—5. We note that also 
the monochromatic flux of Ha at —3000 kms~^ in 
Fig. |9] shows a similar flattening. This velocity is 
dominated by the reverse shock and is subject to 
less systematic errors from the interpolation be- 
tween the lines. 

A similar flattening as seen in the reverse shock 



flux evolution has been seen in other aspects of the 
ring collision. In the optical range the flux evo- 
lution of the narrow lines from the shocked ring 
show a fl attening of most emission lines from day 
~ 7000 (jGroningsson et all [2008ah . These lines 
are formed behind the slow shocks transmitted 
into the ring, mainly as a result of the photoioniza- 
tion by the soft X-rays from these shocks. Their 
flux is therefore a measure of the ejecta - ring in- 
teraction. 



In the X-ray range, iPark et al.l (|201lf ) flnd a 



similar flattening in the 0.5 - 2 keV X-ray flux. The 
hard 3-10 keV flux, however, continues to in- 
crease. Decomposing t he soft flux within the re- 
flected shock model bv iZhekov et al. ( 2010( ). Park 
et al. find a flnd a dramatic decrease of the X-ray 
flux produced by the shock that is transmitted into 
the ring. It is not clear that this is consistent with 
the evolution of the optical lines. The degree of 
this flattening has , however, been challenged by 



Maggi et al 



(|2012f) . The radio emission, corre- 
lated with the hard X-rays, contin ues to increase 
up t o at least 2010 (day 8014) (jZanardo et al 



201fl[ ). 

The ionization of the hydrogen close to the 
reverse shock depends on both the ionizing flux 
and the density. The latter can be estimated 
from models of the ejecta used for light curve 
calculations and spectra d uring the flrst years af- 
ter ex plosion. From the Shigevama fc Nomotol 
()l990l) 14E1 model we can approximate the den- 
sity proflle at 20 years with p{V}_ = 2.0 x 
10l!!(V20 yrs 



-3 V4-«-« g cnr 



Michael et al 



(|2003[ ) flt the lOH model of IWooslevI (|l988h with 



a similar power law, but the density is a factor of 
3.0 higher. The uncertainty in the density is con- 
siderable, both from differences in the explosion 
models and also from the variations in the exact 
position of the reverse shock. We therefore intro- 
duce a factor k to take this into account, which is 
expected to be in the range I ^ k < 3. With an 
H : He abundance of 1 : 0.25 by number we get 
for the envelope density 



i{r) = 71fc [ — 



-8.6 



t 



20 yrs 



5.6 



cm 



-3 



(1) 



where Rs is the radius of the reverse shock. The 
hydrogen density at the reverse shock therefore 
increases rapidly with time as long as the ejecta 
are in the steep part of the density proflle (i.e. 
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V > 4000 kms~^). At lower velocity the density 
gradient flattens considerably, and consequently 
the density at the reverse shock will increase more 
slowly and ultimately decrease as the flat part 
of the density profile is encountered for V < 
3000 — 4000 kms~^. This will, however, only occur 
at > 40 years. This assumes an ejecta structure 
similar to that of the ID models discussed above, 
and instabilities may greatly change this (see be- 
low). 

From 2005 Jan. (day 65 33), which w as the 
epoch of the observations bv ISmith et al. I (|20051 
to 2010 Sept. (day 8619) the 0.5 - 2 keV X-ray 
flux h as increased by a factor of ~ 3.6 ( Park et"al1 
20111 ). From Eq. ([T]) we, however, find that the 



density at the reverse shock during this period has 
increased by a factor of ~ 4.7. The ionization of 
the pre-shock gas should therefore be nearly con- 
stant, or even decreasing in the last epochs. This 
applies especially to the period when the X-ray 
flux levels off. We do therefore not expect a rapid 
decrease of the flux from the reverse shock as was 
predicted on the basis of a steadily increasing X- 
ray flux. 

In the case of a stationary reverse shock, as may 
be the case for the shock in the ring plane, we 
expect the intensity of the Ha line to be 



I. 



in 



"■Hl(^s)Kjccta(R.! 



(2) 



where e = 0.2 — 0.25 is the number of Ha photons 
per ionization between the shocked electrons and 
the free-streaming H I and T 4iecta(Rs) is the eject a 



velocity at the shock (e.g., iMichael et al.l 120031) 



With the density from E q. (U) we find th at /, 



, oc 

t^-^Kijecta OC i''-^ (see e.g.. lHe^ et~alll2006l ) . From 



Fig. [5] we find that the flux of the reverse shock 
has increased by a factor 4.5 from day 5040 
to day 7170. After this it stays nearly constant. 
From the above scaling we would expect it to have 
increased by a factor of ^ 5.1, i.e., close to what 
is observed. 

The scaling above only applies to a stationary 
shock of constant surface area. In reality, the re- 
verse shock is expected to expand (in the observer 
frame) above the ring, where it sweeps up the 
constant, low den sity H II region. In this case 
iHeng et all (|2006[ ) find that the flux should in- 
crease as oc i(2n-9+4s-ns)/(n-s)^ whcrc 71 is the 
power law index of the ejecta density and s that of 



the circumstellar medium. For n = 8.6 and a con- 
stant density medium, we find F ex. The ob- 
served fiux should therefore evolve somewhere be- 
tween these limits, although the ring plane should 
dominate and a flux evolution closer to F oc f^'^ 
is expected, which also gives the best agreement 
with the observations up to ^ 7100 days. 

The nearly constant flux after 7000 days may 
be a result of an early encounter with the higher 
density inner hydrogen region, close to the core. 
The flux evolution will then be determined by 
the area involved in this and the density, F oc 
n{t) po(Vcjocta(-Rs)) t~^, where n{t) is the solid 
angle of the 'high density' encounter and po is the 
density at a reference time tQ. From the HST ob- 
servations in L13 it is indeed found that there is 
especially in the southern part of the ring a bright 
region in Ha, which may come from such a high 
density region of the ejecta. At the last HST ob- 
servation discussed in L13 at 8328 days (22.5 yrs) 
this feature is close to the projected position of 
the reverse shock. Based on its low radial velocity 
it should be close to the plane of the sky, rather 
than the ring plane. Taking a maximum radial 
velocity of 1500 kms~^, the position corresponds 
to an ejecta velocity of ~ 4500 — 5000 kms~^. It 
is therefore likely that the ejecta has a density at 
lower velocities considerably different from the ID 
models. This is a lso confirmed by the simulations 
by Ham mer et al . (2010). An evolution of the re- 
verse shock flux departing from the above simple 
scalings should therefore not be surprising. 

The line proflles combined with the spatial in- 
formation provide important constraints on the 
reverse shock geome try, as was demonstrated by 



Michael et al.l (|2003l ). This applies especially to 
the maximum velocity, which is limited by the 
equatorial ring in the ring plane. 

Our late spectra can give additional informa- 
tion on this from the h igh velocity wings of Ha . 



Sugerman et al.l (|2005f ) find a semi-major axis of 
the ring of 0.82" and an inclination of 43°. This 
gives a radius of the ring equal to6.1xl0^^ cm for 
a distance of 50 kpc. Assuming a radius for the 
reverse shock equal to ~ 80% of the blast wave 
(ring) (|Chevalieii 119821 ) the radius of the reverse 
shock is ^ 4.9 x 10^^ cm, and the ejecta velocity 
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at the shock therefore 



14jccta(i?s) = 7740 



20 yrs 



km s-^ , (3) 



where t is the time since explosion in years. The 
maximum LOS velocity of the ejecta in the ring 
plane is therefore expected to be 



Knax = 5, 660 



t 



20 yrs 



kms~^. (4) 



The fact that we observe velocities up to ^ 
11,000 kms^^ for Ha even in our last 2009 ob- 
servations (Fig. [5]), suggests that the expansion of 
the ejecta must now be anisotropic, with a larger 
extension above the ring and near stagnation in 
the plane of the ring. 

It is interesting to compare our observed Ha ex- 
pansion velocity with the morphology found from 
the modeling of radio ob s ervat ions, as well as Lya 
observations. iNg et all (l2008l) find fro m the ra- 
dio emission that the best fit to the morphology 
is provided by a torus model with a thickness of 
±40°from the equatorial plane, which is similar 
to th e thickness found by I Michael et al. ( 1998b . 
20031 ) from the Lya modeling, ~ ±30°. With an 
inclination of ^ 43° of the ring, this means that 
the reverse shock extends up to nearly the direct 
LOS to the center of the supernova. A maxi- 
mum observed velocity of 11,000 kms~^ there- 
fore also corresponds to a similar radial velocity 
in this direction, and an extension of the ejecta to 
~ 6.9 X 10"(t/20 yrs) cm, or ~ 15% larger than 
the radius of the ring and ~ 40 % larger than the 
reverse shock location in the ring plane. That the 
maximum velocity comes from the part of the re- 
verse shock expanding in the LOS is confirmed by 
the STIS observations in Fig. [TH which shows the 
highest velocity to occur close to the center of the 
ejecta. 

We find little change in either the maximum ve- 
locity or the fiux of the red and blue wings above 
~ 8000 kms^i (Figs. [5]and[ini). The maximum 
velocity is, however, difficult to determine accu- 
rately because of blending with the [O I] lines from 
the ejecta on the blue side and narrow lines on the 
red. The STIS observations in Fig. [191 however, 
show that there does seem to be a decrease of the 
reverse shock velocity in Ha close to the LOS of 
the center of the ring. Although contaminated by 



the [O I] emission from the ejecta especially in 
2004, we find in the 2004 spectrum a maximum 
blue velocity of -12,200 kms^^, while in 2010 
this velocity has decreased to —9100 knis^^. The 
corresponding values on the red side are -1-8300 
kms~^, and -1-7400 kms^^. The velocities on the 
red side are lower than the maximum measured 
with UVES and F0RS2, which are likely to be a 
result of the limited S /N of the STIS spectra. 

The steepening of the line profile with time is 
roughly consistent with the maximum LOS veloc- 
ity of the ejecta in the plane of the equatorial 
ring (see Sect. 14. 1[) . This is the point of maxi- 
mum interaction, and it is therefore not surpris- 
ing to find that the most dramatic change oc- 
curs aro und this veloc i ty, as was concluded al- 
ready by I Michael et al. ( 2003 ) . If the interaction 
takes place in a thin cylinder concentric to the ring 
with small extent in the polar direction, and with 
roughly equal strength along the rim, it will re- 
sult in a double peaked line profile between ±Vmax 
(Eq. |4]). As the thickness of the shell increases in 
the polar directi on the line approaches a flat, boxy 
profile (see Michael et al.l i2003. for a discussion). 
The increasing flux of the central peak of the line 
is therefore partly a result of the increase of this 
component. As Figs. [8] and |9] show, there is, how- 
ever, also a substantial increase of the flux from 
the inner ejecta component on top of this. 

The fact that there is emission above ^ 8000 
kms~^ may require the ejecta to interact with 
a medium of comparable density to that in the 
plane. At 20 years and 11,000 kms~^ we get from 



Eq. ([TJ a number density 



3 cm 



With 



this density profile we find a ratio between the 
ejecta density in the ring plane, where V^j w 7740 
kms~^, and the density at 11,000 kms^^ equal 
to ^ 20. With an H II de nsity of 100 cm'^ 
(jChevalier fc Dwarkada£lll995l) an estimate based 
on 14 f« {Pei/pmiY^'^Vej gives a shock velocity of 
^ 1900 kms~^ into the H II region above the ring 
plane. 

According to lMattila et all (|2010[) , 100 cm"^ H 
II region gas can, however, not fill the full volume 
from the ring plane up to our LOS, as this would 
violate observational constraints on the narrow 
line emission. It is therefore likely that the H II 
region is less dense along our LOS than 100 cm~'^, 
which results in a higher shock v elocity than 1900 
km/s. This is also indicated by iNg et al.l (|2008h 
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who find from the radio a mean expansion speed 
of 4000 ±400 kms-i over the period 1992 - 2008. 
However, X-ray observations by Chandra indicate 
a consi derably lower exp ansion speed of ^ 1625 
kms-i (|Park et al.ll201l]) . 



The high velocity emission from the ejecta in 
the LOS could also be a result of the X-ray ion- 
ization from the ring. If a large fraction of the 
X-ray flux is below ^0.1 keV these X-rays may 
be absorbed in the outer ejecta and at high alti- 
tudes from the ring (see Fig. [22] below). They 
can there provide the necessary energy for the Ha 
emission. The fact that there is little change in 
the flux in the high velocity wing of Ha, however, 
argues against this explanation. 

Adopting a similarity solution we expect Vcj oc 
giving a decrease in the ejecta velocity by ~ 
20% over 10 years. There is reason to doubt that 
the similarity solution applies over the small time 
range since the impact on the H II region. With a 
maximum velocity of 11,000 kms~^ this may be in 
conflict with the observations (Fig. [5|) . As pointed 
out in Sect. 13.11 the maximum velocity could be 
higher and a decrease would then be 'hidden' in 
the many lines from the ring collision. 

4.2. Ejecta component 

Based on photometry of the ejecta with the 
WFPC2, ACS and WFPC3 on HST, Lll find that 
the flux has increased by a factor of 2 - 3 in the R- 
and B-like bands from day 5000 to day 8000. The 
R-band is dominated by Ha, wh ile the B-band is 



a mi x of H/3 and Fe I-II lines (|Jerkstrand et al 
201 ih . 



Our light curve of Ha (Figs. [8] and [9]) shows an 
increase by a factor of 4 — 6, while the [Ca II] A A 
7292, 7324 lines (Fig. ^ increased by a factor of 
5 — 8 over the same period. Given that the HST 
observations cover a wide spectral range for each 
filter and that there is a considerable uncertainty 
in the line fluxes because of blending with the 
narrow lines, this is consistent with the increase 
found from the HST broad band photometry. This 
therefore conflrms spectroscopically what was con- 
cluded based mainly on imaging in Lll. In addi- 
tion, we see that other lines, in particular the Mg 
II AA 9218, 9244 lines (see Fig. |4|), show a simflar, 
or even larger increase. 

Later than - 7000 days the Ha, [Ca II] and Mg 



II fluxes level off, similarly to the reverse shock 
lines and the X-rays, discussed in the previous 
section. As we discuss in next section (see also 
Lll), there is a close connection between the X- 
rays from the ring collision and the ejecta flux. It 
is therefore expected that the ejecta flux follows 
the flattening seen in X-rays. 

As for the line proflles from the ejecta. Ha is 
severely contaminated by the narrow lines, as well 
as the reverse shock emission. The [Ca II] lines 
are less affected by these and also lack a reverse 
shock component (Fig. [T8|) . The line shape on 
the blue side is well deflned and approximately 
triangular, similar to that of the Ha line in the 
early STIS observation in Fig. [lH The blue side 
reaches ^ 4500 kms~^, and does not change ap- 
preciably from 2002 to 2012 (days 5738 - 9019). 
This velocity is similar to Ha, as measured from 
the 1999 STIS G750M spectrum in Fig. [H Both 
the similar line profiles and the similar maximum 
velocities indicate that these lines arise from the 
same region. 

L13 show that the changing morphology results 
from the same cause as the increase observed in 
the light curve. They find that the Ha emission 
from the inner ejecta has changed from a centrally 
dominated, elliptical profile before year 2000 to 
one dominated by a ring-like structure at later 
stages. Given this complex morphology it is some- 
what surprising that the line profiles from the core 
component change relatively little, except in flux 
(see Figs. [5] and [T2|) . As noted above, a reason 
for this may be that most of the change occurs for 
LOS velocities less than ^ 1500 kms~^, a spectral 
region which as discussed earlier is blocked by the 
narrow lines. 

There is also little evolution of the [Si I]/[Fe II] 
line from the VLT/SINFONI observations (L13). 
This emission is, however, centrally peaked and is 
likely to come from more central, processed mate- 
rial, protected from most of the X-rays, and does 
not display any major change in morphology be- 
tween 2005 and 2010. 

4.3. The energy deposition by the X-rays 
from the ring collision 

Lll show that up to day 5000 the light curves 
in the R and B-bands are compatible with the ra- 
dioactive decay of ''^Ti. They concluded that to 
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explain the increasing ejecta flux an additional en- 
ergy source is needed. The most likely such source 
is the soft X-rays emitted by shocks resulting from 
the ejecta - ring interaction. We now discuss the 
details of the deposition of these and their conver- 
sion into optical/UV radiation. 

The observed X-ray luminosity above 0.5 keV 
is ~ 4 X lO^^ergs" ^ at the last observation at 8619 



davs bv iPark et al 



Based on the spectra, 
Zhekov et al.l (|201Cll ) has modeled the X-ray emis- 
sion with three components, one fast blast wave, 
one soft component from the transmitted shocks 
into the dense clumps in the ring and one from 
the reflected shock component. The soft radiation 
from the transmitted shocks into the ring, with 
kTe ^ 0.35 keV, dominates the total flux. 

After ~ 10 years the hydrogen envelope was 
mostly transparent to X-rays with energy > 1 
keV (Lll). The inner core region, containing 
both mixed-in hydrogen and synthesized heavy el- 
ements, which increase the opacity, is, however, 
still opaque up to ~ 8 keV. The X-rays will there 
be thermalized into optical, IR and UV radiation. 
Only after > 30 years does the whole core become 
transparent. 

The exact X-ray luminosity and spectrum is un- 
certain. Observations only probe energies down to 
~ 0.5 keV. For low shock velocities most of the X- 
ray flux is below the lo w energy interstellar X-ray 

(I2OO8H) 



cut-off at ~ 0.5 keV. In Groningsson et al 



it is shown that the shocked lines have a FWHM 
~ 200 kms~^. Even if there could be an incli- 
nation effect, which could increase this velocity 
by a factor up to ^ 1.4, the typical shock veloc- 
ities giving rise to these lines are < 300 kms~^, 
correspondi ng to shock temperatures < 0.1 keV. 
In addition, iGroningsson et al.l (|2008al) show that 
shocks into the densest parts of the ring are radia- 
tive up to ~ 500 kms~^, or ~ 0.25 keV. Therefore, 
a large fraction of the X-rays/EUV flux may well 
be below the observed ISM absorption threshold 
at ~ 0.5 keV. 

The X-rays give rise to secondary electrons 
by photoionization in the same way as the 
amma-rays and po sitrons from rad ioactivity 
Shull fc van Steenberg 1985; Xu fc McC rav 199^; 



Kozma fc Fransson 19921 ). These in turn will 



lose their energy in non-thermal excitations, ion- 
izations and Coulomb heating of the free ther- 
mal electrons. The resulting optical/UV spec- 



trum from the X-ray input will therefore not be 
very different from that of the radioactive de- 
cay, as long as the ionization of the gas is low. 
For a low ionization each hydrogen ionization re- 
quires ^ 30 eV at th e relevant level of ionization 
(|Xu fc McCravl Il99l[ ) . Once the ionization in- 
creases above ~ 10~^ the efhciency of excitation 
and ionization decreases, while that of the heating 
increases. The cooling, balancing the heating, is 
close to the core region at these epochs mainly 
done by thermal, collisional excitations of mid- 
and far-IR fine structure lines and molecules, and 
possibly cool dust. This emission will therefore 
be difficult to observe directly, although far-IR 
and mm observations may help. In the envelope 
adiabatic cooling dominates. 

Compared to the radioactive excitation some 
important differences arise due to the fact that 
the X-ray illumination is from the outside, while 
that of the gamma-rays and positrons is from the 
central regions. The positrons from the ^^Ti de- 
cay, which dominate the radioactive energy in- 
put at these epochs, wiU be deposited locally, 
mainly in the Fe-rich gas, unless the magnetic field 
in the ejecta is extrem ely low (see discussion in 



Jerkstrand et al.l 120111 ). The X-rays will on the 



other hand penetrate to different depths, depend- 
ing on their energy. 

The photoelectric energy deposition is sensitive 
to the density and abundance distributions in the 
ejecta. To calculate this we use the eje cta model, 
14 E1 from Shigev ama fc Nomotol ( 1990l) . mixed as 



Blinnikov et al.. (.2000. ). This is a spherically 



symmetric model, where the chemical mixing has 
been modified to reproduce the light curve. In re- 
ality, Rayleigh- Taylor instabilities will mix the dif- 
ferent nuclear burning zones in velocity, resulting 
in both l arge abundance and density inhomogen- 
ites (e.g.. Hammer et al. 2010h . Nevertheless, this 
model should show the main qualitative features 
of the X-ray deposition. 

Using the same code as was used in Lll for 
calculating the total energy deposition of the X- 
rays, we can calculate also the distribution of the 
energy deposition in radius. We assume here that 
the X-ray sources are located in the plane of the 
equatorial ring at 6.1 x 10^^ cm and the ejecta 
extend to the reverse shock at Rs = 4.9 x 10^^ cm. 
In addition, we assume spherical geometry for the 
ejecta, which is certainly an oversimplification. 
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In Fig. [21] we show the fractional energy de- 
position per cm integrated over the polar angle, 
df /dr, for energies between 0.1 - 5 keV at 20 and 
30 years after explosion. The total energy depo- 
sition in a shell with radius r and thickness dr 
is then dL{E) = Lx-ray(£^) df{E)/dr dr. As we 
show below, this is, however, only a rough approx- 
imation and the real deposition function is two di- 
mensional or even three dimensional. 

We here see that for the hard X-rays with > 1 
keV most of the energy is deposited at the outer 
boundary of the core, while the softer X-rays are 
deposited in the hydrogen envelope. At 30 years 
the soft X-rays are deposited in a more narrow 
velocity interval than at 20 years, caused by the 
increasing ejecta density close to the reverse shock 
(Eq. [T|). The harder X-rays, however, are de- 
posited deeper in the core, where the flatter den- 
sity distribution does not have the same effect. 
The fact that the core is very heterogeneous in 
abundance, with a large fra ction of the v olume 
occupied by the iron bubble (jLi et al.lll993l ). hav- 
ing a large X-ray opacity, means that most of the 
hard X-rays will be absorbed efficiently here. 

Because the interaction of the ejecta and ring 
is mainly taking place in a thin region around 
the equatorial ring, the X-ray illumination of the 
ejecta will be highly non-uniform in the polar di- 
rection. Most of the energy deposition will oc- 
cur close to the equatorial plane, but especially 
at higher energies a large fraction of the ejecta 
will be transparent and the deposition more dis- 
tributed also above the plane. To calculate the de- 
position in 3-D we assume that the X-rays are in- 
jected in a thin region in the equatorial plane, and 
azimuthally uniform. This is a reasonable approx- 
imation for observations obtained after ^ 2004, 
but for earlier epochs the distribution of the indi- 
vidual hotspots in the azimuthal direction should 
be taken into account. In Fig. [22] we show the de- 
position at different energies in the poloidal plane. 

Again, we note the larger penetration of the 
hard X-rays. More interesting is, however, that 
the X-rays are deposited spatially fairly locally. At 
low energies this is close to the ring and also close 
to the equatorial plane, while the higher energies 
are deposited in a more extended region closer to 
the core. We also note the higher concentration 
of the energy deposition to the ring at 30 years, 
compared to 20 years. 
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Fig. 21. — Energy deposition per cm for the 
mixed 14E1 model for different photon energies 
at 20 years and 30 years after explosion. We 
also show t he density distribution from the 14E1 
model from Blinnikov et al. I (I2OOOI) scaled to these 
epochs. Note the sharp density gradient in the 
core at r ^ 1.0 X 10"'^'' cm and the power law gra- 
dient for r > 1.5 X 10^^ cm (both at 20 years), 
corresponding to ^ 4000 kms^^. The velocity at 
the ring is 10,500 kms^^ at 20 years and 6900 
kms~^ at 30 years. 



As mentioned in Lll, the restricted radial range 
of the energy deposition may affect the morphol- 
ogy of the HST images, which show a ringlike 
structure in the ejecta. This may be explained as 
a result of the limited penetration of the X-rays 
from the shock, and therefore reflects the energy 
deposition in the ejecta rather than the true den- 
sity distribution, as was seen in Fig. [5TJ This is 
discussed in more detail in L13. 

Figure [21] also shows that if most of the X-ray 
luminosity is below ~ 0.5 keV, as may well be the 
case, the outer parts of the ejecta may absorb a 
larger fraction of the energy, while the core only 
receives a minor fraction of the X-ray input. A fur- 
ther brightening of the outer parts can therefore 
be expected. As has been ment i oned earlier, re- 
cent observations bv IPark et aD (|201ll ). however, 
show a flattening of the soft (0.5 - 2 keV) X-ray 
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Fig. 22. — Energy deposition for the same model 
as in Fig. [21] but now shown perpendicular to the 
equatorial plane at 20 years (upper panels) and 
30 years (lower panels) after explosion. The x- 
and y-scales are in lO^*" cm. The red circle at 
6.1 X 10^^ cm marks the position of the source of 
X-rays (i.e., the ring), while the blue circle with 
radius 4.9 x 10^'' cm marks the location of the re- 
verse shock, assuming spherical symmetry. The 
different contours and colors are distributed ac- 
cording to the logarithm of the energy deposition 
per volume, shown in the color bar at the bottom 
of each panel. 



flux (see however iMaggi et al. ( 2012 1) ) , while the 
hard (3 - 10 keV) flux continues to increase. 

In the high velocity hydrogen envelope the den- 
sity is low, and the recombination timescale long. 
For a density profile similar to that in Eq. ([Ij the 
density at the reverse shock is 70 — 200 cm~'^ at 
an age of 20 years. In the a bsence of X-ray in- 
put w e find using the code in Fransson fc Kozinal 
(|2002l ) that the temperature is here expected to 
be ~ 10 — 20 K . E xtrapola t ing th e Case B recom- 
bination rate from iMartinI ( 19881 ) gives a recom- 
bination rate ~ 2 x 10~^^ cm^s~^ at 20 K and a 



22 kr 



8.6 



20 yrs 



-5.6 



years. 



(5) 



where Xe is the ionization fraction. Equation ([5]) 
shows that the recombination time decreases fast 
both with time and with decreasing radius. As 
long as <C 1 it will, however, be long at least 
close to the shock. As the ionization increases re- 
combination becomes more and more important. 

From the same calculation as above, again in- 
cluding only the radioactive input, we find that 
the degree of ionization in the envelope is nearly 
constant with time at ~ 7 x 10~^, as the 
ionization is frozen-in. Because of the X-ray il- 
lumination the ionization may, however, increase 
considerably. Because of this and the large den- 
sity gradient, the recombination timescale may be- 
come short only a short distance inside the reverse 
shock. A stationary situation may then be at hand 
where absorption and emission balance. 

For the regions where the recombination timescale 
is long, the ionization of the envelope is deter- 
mined by the number density of H atoms in the 
part of the envelope where the X-rays are de- 
posited and the time integrated ionizing luminos- 
ity, rather than the instantaneous luminosity. Fur- 
ther in, once the ionization has increased enough 
there may be a balance between ionizations and 
recombination. 



Int egrating the X-ray luminosity from lPark et al.l 
POIII ) over time gives a total X-ray energy 
Ex-vay ~ 1-6 x lO"*-* ergs in the 0.5-2 keV band. 
Two thirds of this is emitted in the period 19 - 22 
years after explosion. On top of this there is a mi- 
nor contribution from hard X-rays and a possibly 
dominant contribution below 0.5 keV. Roughly 
half of the X-ray energy will be emitted inwards 
to the ejecta. As shown above, most of the X-rays 
and EUV photons below ~ 0.3 keV are absorbed 
by the envelope. The luminosity below this energy 
is therefore most important for the ionization of 
the envelope, but also the most uncertain obser- 
vationally. 

To take into account the fact that we do not 
know the level of the unobserved soft X-ray flux, 
we assume that the total time integrated luminos- 
ity emitted in a narrow energy interval in the EUV 
is E'euv = / X 10"*^ ergs, where / is a scaling 
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factor reflecting the uncertainty in the EUV flux. 
To get an approximate estimate of the ionization 
in the ejecta we further assume that recombina- 
tions are slow and that each volume is absorbing 
a fraction of E^vv given in Fig. [251 We also ig- 
nore the decreasing efficiency of ionization as the 
ionization increases. Finally, for simplicity we ig- 
nore the expansion, which is justified by the fact 
that most of the X-rays are emitted during a pe- 
riod short compared to the age of the supernova. 
These assumptions affect the quantitative results, 
but hardly the qualitative. With these assump- 
tions the ionization of a given volume depends on 
the deposited energy and is inversely proportional 
to the density. 

In Fig. [23] we show the resulting distribution 
for two values of the efficiency of ionization, / = 1 
and / = 0.1, which may bracket the range. Not 
surprising, the distribution reflects the energy in- 
put in Fig. [221 More interesting is that close to the 
shock the ionization is expected to be fairly high, 
in particular for the case where most of the ion- 
izing flux is below ~ 0.3 keV. Comparing the 0.1 
keV and 0.35 keV cases, we note the higher, but 
more concentrated ionization in the former case 
for the same total energy. 
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Fig. 23. — Ionization for the same model as in Fig. 
[22] at 20 years after explosion. The different con- 
tours and colors are distributed according to the 
logarithm of the fractional ionization, Xe, shown in 
the color bar at the bottom of each panel. The two 
sets give the ionization for two different assump- 
tions of the ionizing flux. The x- and y-scales are 
in 10 
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comes from direct excitation of these. The frac- 
tion of the deposited energy going into excitation 
of levels with n > 2 is ~ 12% for Xe = 10"^ 
()Xu fc McCravlll99ll) . Of this ~ 10% is emitted 
as Ha emission. Assuming that we underestimate 
the UVES luminosity by a factor of ~ 2, as was 
found from the 2004 STIS observation (Sect. K^ . 
our observation at ~ 8000 days (Fig. [8]) result 
in L(Ha) « 3.2 x lO^" ergs'^ To explain this 
with direct excitation only, an X-ray luminosity 
of ~ 102L(Ha) « 3.2 x lO^^ ergs"! is therefore 
needed. In the regions where the recombination 
time may become short compared to the age also 
the non-thermal ionizations and subsequent re- 
combinations will contribute to the Ha luminosity, 
increasing this by a factor of ~ 3, compared to the 
direct excitations. The observed 0.5-2 keV X-ray 
luminosity at 800 days was ~ 3.6 x lO'^^ergs"^ 



If the recombination time close to the reverse 
shock is long most of the emission in the H I lines 



([Park et al.ll201l[ ). An additional 14% was in 
the hard 3-10 keV range. There does therefore 
seem to be enough X-ray luminosity to explain the 
observed Ha luminosity. 

The conclusion of this discussion is that the ion- 
ization in both the envelope as a whole and the 
core are likely to remain low for at least the next 
ten years. A dramatic change of the spectrum is 
therefore not expected, although the line ratios 
may slowly change. 

The estimate bv lSmith et al. ( 2005 ^ of the ion- 
ization assumes that most of the X-rays are ab- 
sorbed close to the reverse shock. The H I density 
in front of the reverse shock is at 20 years ~ 
90 cm^'^. The mean free path close to the Lyman 
limit is then - 1.6x 10i''^(£'/13.6 eYf{n/100 cui-^)' 
cm, so most of the ionizing photons are indeed ab- 
sorbed close to the reverese shock as long as the 
photon energies are not too high. As Figs. [2T1 
and [22] show, this is in general the case for photon 
energies < 0.1 keV, but not for higher energies. 
The fact that we still observe both Ha and Lya 
shows that the X-rays are not yet able to pre- 
ionize the gas completely. As we see from Fig. l23] 
the ionization decreases substantially above the 
ring plane, and at least a fraction of the Ha may 
originate here. 

Except for the Mg II AA 9218, 9244 lines, there 
are no major changes in the relative fluxes of the 
lines from the core (Fig. l4]). The Mg II lines, 
however, increase substantially compared to e.g., 
the [Ca II] lines. If these lines are pumped by Lya 
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there are several factors which influence the outgo- 
ing flux. The most straightforward is the increas- 
ing Lya flux, which should be reflected directly 
in the Mg II lines. However, also the fractional 
population of Mg II, which determines the optical 
depth of the Mg II lines, and therefore the pump- 
ing efflciency, may increase because of the ionizing 
effects of the X-rays. Finally, the fraction of the 
envelope which is in velocity resonance with the 
pumping transition may change with time. 

4.4. Dust 

Strong evidence for early dust formation in the 
ejecta w as found from bo th blueshifts in the line 
profiles (|Lucv et aLlll989l) and a thermal IR excess 



( Wooden et al.Ml993D . The last far-IR observa- 



tions before the ring i n teract ion are from day 1731, 
where Bouchet et al. ( 19961 ) find a dust tempera- 
ture of ~155 K. There is also from the line profiles 
of Mg I] A 4571 some evidence for dust from the 
HST observations at 8 ye ars with an effective op - 
tical depth close to unity (jJerkstrand et al.luOllI) . 
This is, however, likely to be in the form of very 
dense clumps with high optical depth. The effec- 
tive optical depth is therefore only a measure of 
the covering factor of the dust clumps. 



More recentlv lMatsuura et al.l (|2011l ) have from 
Herschel far-IR observations found evidence for a 
large dust mass, estimated to 0.4 — 0.7 Mq, but 
could be up to 2.4 Mq. The mass estimate is, how- 
ever, highly dependent on the dust temperature 
and the mass could be several orders of magnitude 
lower if the te mperature is higher than adopted by 
Matsuura et al. (2011 ), as discussed f rom APEX 
observations by Laki cevic et al.l (|2012r) . 

As was discussed in Sect. 13.21 the 1999 Ha line 
profile may indicate a ^ 15% deficit on the red side 
below ^ 2000 kms^^ (Fig. ITS]) . This only repre- 
sents weak evidence for dust in the ejecta. The 
2D spectra do, however, show a clear asymmetry 
with most emission on the blue side. While the ob- 
served asymmetrical emission in Ha is a possible 
explanation the prevalence of blue emission is in- 
dicative of dust obscuration. Because most of Ha 
arises outside the metal core, where the dust pre- 
sumably is located, it is, however, less affected by 
dust in the metal core than e.g, the [O I] emission. 
There are also some indications of asymmetries 
from the [Ca II] line profile (Fig. [T51) . consistent 
with that in Ha. As already mentioned, this line. 



however, originates in the same region as Ha. 

The most direct evidence from the line profiles 
at early epochs came f rom the evolution of the [O 
I] AA 6300, 6364 lines (|Lucv et al.lll989t ). Lucy et 
al. found a blueshift of the peak by ~ 600 kms~^ 
at the time of the peak velocity shift, ^ 530 days. 
Figure [11] represents the last observation of this 
line before the spectrum becomes affected by the 
X-rays. We here note that the peak of the [O I] A 
6300 line also in this late observation is blueshifted 
by 1000 kms~^. The resolution of this spectrum 
is limited, but the blueshif t is con sistent with that 
determined by Lucv et al. ( 19891 ). 

L13 argue that the 'hole' seen in the optical 
HST images is the result of the external X-ray 
illumination rather than dust obscuration. We 
can, however, not exclude that this region con- 
tains some dust. After all, this metal rich re- 
gion is the most likely dust forming region. The 
VLT/SINFONI IFU spectra from 2005 and 2010 
also show a smooth decrease of the red wing of the 
[Si I]/[Fe II] A 1.644 /xm line at velocities < 3000 
kms""'^ (L13). This is roughly what is expected 
from dust internal to the line emitting region. 

Summarizing this discussion, we find indica- 
tions of dust from the 1999 [O I] A 6300 line (and 
also the [Si I]/[Fe II]) profile, but also caution that 
the asymmetries in the core may also give this ef- 
fect. From the high resolution observation of the 
Ha line there is, however, little indication of dust, 
although this line should be less affected by this. 

5. Summary and Conclusions 

We have here discussed the spectral evolution 
of the ejecta and reverse shock emission from 2000 
to 2012 (days ~ 4400 - 9100). Together with the 
HST imaging this represents a unique data set. 
Our main conclusions are summarised below: 

• Both the Ha, [Ca II] and Mg II lines from the 
inner ejecta have increased by a factor of 4-6 
from 2000 to 2012. This confirms the broad 
band increase found from the HST observa- 
tions by Lll. There is a flattening of the 
light curve later than ^ 7000 days, possibly 
correlated with a similar flattening of the X- 
ray flux and the flux from the narrow lines 
from the shocked ring. 

• The reverse shock flux behaves in a similar 
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way. The leveling off of the flux is probably 
not a result of an increasing ionization of the 
pre-shock hydrogen, but may be a result of 
an arrival of protrusions of the inner, dense 
core region at the reverse shock, as seen from 
the HST imaging. 

The reverse shock emission extends to veloc- 
ities > 11,000 kms~^. This is larger than 
can be contained within the equatorial ring 
and indicates an anisotropic expansion of the 
ejecta above and below the ring plane. 

The Balmer and metal lines have an inner 
ejecta component with a width of ~ 4500 

kms""'^. 

From the last radioactively dominated STIS 
spectrum in 1999 (day 4381 - 4387) we find 
similar velocities for the Ha, Na I and [Ca II] 
lines, consistent with an origin in the hydro- 
gen rich gas. The [O I] has a lower velocity, 
mainly coming from the core. 

There is from Ha and [O I] A 6300 evidence 
of a blueshift similar to what was seen as a 
dust indicator at early epochs. 

We identify a line feature at ^ 9220 A with 
an Mg II line pumped by Lya fluorescence. 
The Mg II omission is mainly from primor- 
dial magnesium in the H and He rich ejecta. 
This line has been increasing faster than the 
other lines from the inner ejecta. 

The deposition of the X-rays from the shock 
region affects mainly the hydrogen envelope 
and outer core region of the ejecta. Only 
X-rays with energy > 5 keV reach the inner 
core region. 

The changes in the morphology of the ejecta 
during the last decade is probably mainly a 
result of the increasing ionization and heat- 
ing of the ejecta by the X-rays. The steep 
density gradient in the hydrogen envelope 
may lead to a more concentrated ionization 
to the reverse shock with time. 

The excitation by the X-rays is similar to the 
gamma-ray and positron deposition as long 
as the ionization is low. Most of the heat- 
ing is balanced by far-IR fine structure and 



molecular lines, while the optical and near- 
IR lines are powered by non-thermal excita- 
tion and ionization. 

Future observations are needed to confirm the 
evolution and follow the increasing ionization of 
the ejecta by the external X-rays. 
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